Liquid ^He near the superfluid transition in the presence of a heat current and gravity 
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The effects of a heat current and gravity in liquid ''He near the superfluid transition are investi- 
gated for temperatures above and below Tx. We present a renormalization-group calculation based 
on model F for the Green's function in a self-consistent approximation which in quantum many- 
particle theory is known as the Hartree approximation. The approach can handle the average order 
parameter (i/>) = above and below T\ and includes effects of vortices. We calculate the thermal 
conductivity At (AT, Q) and the specific heat C{AT, Q) for all temperature differences AT = T — T\ 
and heat currents Q in the critical regime. Furthermore, we calculate the temperature profile T{z). 
Below Ta we find a second correlation length ~ Q^^{T\ — T)^'^ which describes the dephasing of 
the order-parameter field due to vortices. We find dissipation and mutual friction of the superfiuid- 
normal fiuid counterfiow and calculate the Gorter-Mellink coefficient A. We compare our theoretical 
results with recent experiments. 
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I. INTRODUCTION 

Gravity and a heat current Q are two sources which in- 
fluence the superfluid transition of liquid ^He at Ta « 2 K 
and cause inhomogeneities in the system. On earth 
gravity implies a pressure variation P — P{z) where z 
is the altitude coordinate. Since the superfluid tran- 
sition temperature T\ — Tx{P) is pressure dependent, 
Tx(z) = rx(P(z))|-|depends on the altitude coordinate z 
with the gradienti dTx/dz = -hi. 273 ^K/cm. On the 
other hand a nonzero heat current Q drives the system 
away from equilibrium. A temperature gradient VT is 
created which implies that the temperature T is space 
dependent. We assume that the heat current Q is ho- 
mogeneous and flows vertically (parallel to the z axis) so 
that the temperature T(z) depends on the z coordinate 
only. 

The local properties of the system are determined 
by three parameters, the local temperature difference 
AT{z) = T{z) - Tx{z), the heat current Q, which is 
related to VT, and gravity g, which is related to VT\. 
The point (AT, Q, g) = (0, 0, 0) is a critical point related 
to the superfluid transition. This means that in thermal 
equilibrium {Q = 0) and in microgravity {g = 0) the 
system shows a second-order phase transition at T = T^ 
from the normal-fluid to the superfluid state. Usually, 
gravity is negligible except for very small heat currents 
Q and for very small AT, i.e. very close to Tx- Since on 
earth the gravity acceleration g = 9.81 m/s^ is a flxed 
quantity, in most cases the AT-Q plane is considered as 
the phase diagram. 

Liquid ^He close to Tx in the presence ofp.a^ heat 
current Q has, been investigated theoreticallytra and 
experimentallyQ~l3. In the AT-Q plane a line of critical 
temperatures is found which separates superfluid from 
normal-fluid helium. A nonzero heat current Q implies 
that the superfluid transition temperature Tx is shifted to 



lower temperatures by ATx{Q). Foiiiimall heat currents 
Q in the critical regime the theoryou predicts the shift 

ATa(Q) Q="j}#^the exponent x = l/2z/ = 0.745. In 

the experimentsErEliliJ a depression ATc{Q) of the super- 
fluid transition has been observed which agrees qualita- 
tively with the theory, but not quantitatively. For small 
Q in the critical regime the shift ATc{Q) ^ — Q^-^^ has 
been founcO. While the exponents do not agree, the ex- 
perimentally observed shift ATc{Q) is larger than the 
theoretically calculated ATx{Q)- 

While for AT > ATx{Q) the hehum is normal fluid and 
for AT < ATc(jQ) it is superfluid, in a recent experiment 
Liu and AhlersO found a new dissipative region for tem- 
peratures AT in the interval ATc(Q) < AT < ATa(Q). 
This observation indicates that at a flnite heat current 
Q the transition from normal fluid to superfluid helium 
may possibly happen in two steps with two transition 
temperatures ATx{Q) and ATc{Q) (relative to the equi- 
librium transition temperature Tx). While the upper 
ATx{Q) may be identifled by the theoretical prediction, 
the lower ATc(Q) agrees with the shift of Ref. In a sim- 
ilar experiment performed by Murphy and Meyeru also 
two transition temperatures were found. While the val- 
ues of ATx{Q) and ATc{Q) contain errors, the difference 
ATx{Q) — ATc{Q) is quite well reproduced by the latter 
experiment. 

Heat-transport phenomena in liquid '*He close to Tx 
are described by model F of Halperin, Hohenberg, and 
Siggialla which is a model for the critical and hydrody- 
namic slow variables including fluctuations. Most the- 
oretical investigations are based on this model. In the 
normal-fluid region for temperatures T above Tx the 
heat is transported diffusively driven by the tempera- 
ture gradient VT. In linear response the heat current is 
Q = — AtVT where At is the thermal conductivity. For 
inflnitesimal Q and zero gravity the thermal conductiv- 
ity At has been calculated within model F in two-loop 
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ordeiEi. Critical fluctuations, which are taken into ac- 
count by the renormalization-group (RG) theory, imply 
a strong enhancement of At close to Tx. For infinites- 
imal Q and zero gravity the thermal conductivity At 
diverges in the limit T — > T\, or more precisely in the 
limit {AT,Q,g) (0,0,0). The RG theory has been 
extendedB to calculate At for nonzero heat currents Q 
but without gravity. It turns out that for finite Q close 
to Tx the heat transport becomes nonlinear which means 
that At becomes Q-dependent. It has been shownB that 
At remains finite in this case even for T = T\. 

On the other hand, in the superfiuid region where 
AT{z) < ATc{Q) the heat is transported convectiveJM 
nearly without friction according to the two-fluid modelcj 
by the superfluid-normal fluid counterflow. In this case 
the temperature gradient VT is nearly zero indicating a 
nearly infinite thermal conductivity At. Mutual friction 
between the superfiuid and the normal-fiuid component 
and dissipation of the heat current occur only by creation 
of vortices, which however is a small effect. Nevertheless, 
mutual friction in the superfiuid state has been measured 
as early as 1949 by Gorter and Mellinktj. For the mutual 
friction force the ansatz / — Ap^psivs — Wn)^ was madcEj 
with a temperature dependent coefficient A, the so called 
Gorter-Mellink coefficieat. This ansatz is related to a tur- 
bulent superfiuid flowta and implies a temperature gra- 
dient VT ~ —Q^, so that the thermal conductivity is 
At ^ Q^^- More recently, the temperature gradient VT 
due to mutual friction ip-|Superfluid '*He was measured 
directly by Baddar et alO in the critical regime close to 
Tx. This experiment confirms the ansatz by Gorter and 
Mellink qualitatively but with a slightly different expo- 
nent in the Q dependence or — '^n dependence. 

In the intermediate region a crossover between the two 
heat-transport mechanisms happens. Since the system 
is spatially inhomogeneous this crossover happens also 
spatially and implies an interface between superfliiid and 
normal-fluid helium located at a certain zq. Onukio inves- 
tigated this interface by solving the model-F equations in 
mean-field approximation, where critical fiuctuations are 
taken into account by scaling theory. He calculated the 
temperature profile T{z) and the order-parameter profile 
ipiz). While in the normal-fluid region the temperature 
profile has a finite gradient related to a finite thermal 
conductivity At, in the superfiuid region the tempera- 
ture profile is absolutely flat and the gradient is zero, so 
that mutual friction is not iiacjljuded and At is infinite. 

The previous RG theorieada were constructed as per- 
turbation theories starting with mean-field solutions of 
the jpodel F equations. While in the normal-fluid region 
it iaJ (ip) — 0, in the superfluid region a plane wave order 
parameter (ip) = rj e**^^ was assumeda. Consequently, the 
temperature profile was found to be fiat in the superfiuid 
region so that mutual friction and dissipation by vortex 
creation are not included in the previous RG theories. 
In this paper we extend the RG theory in two respects. 
First, we include gravity. Secondly, we include mutual 
friction and dissipation in the superfiuid region. 



Strong fiuctuations of the phase of the order-parameter 
field ip can imply that the average order parameter is 
{ip) = even below T\. This fact is well known for 
systems of finite size. Here for nonzero Q the vortices 
cause sufficiently strong phase fluctuations so that (ip) = 
not only above but also below Tx- For this reason, in 
the noncquilibrium state with a nonzero heat current Q 
we may not start a perturbation theory with a nonzero 
mean- field order parameter which implies a nonzero (rp) . 
Rather we need a new approach which can handle (ip) — 
above and below Tx- In this paper we present a self- 
consistent approximation for the Green's function which 
can do this and which in quantum mamtparticle theory 
is known as the Hartree approximationES. 

The paper is organized as follows. In Sec. || we 
describe the model and the necessary field theoretic 
tools, i.e. the Feynman rules. In Sec. |lll| the idea of 



i.e. the Feynman rules. In Sec. 
our approach and the calculations of the Green's func- 
tion and the effective parameters are presented. The 
renormalization-group theory is applied in Sec. ^ to in- 
clude the critical fluctuations. The thermal conductivity 
and the temperature profile are calculated in Sees. [v|and 
VI, while in Sees. VII and VIH the correlation lengths, 
the entropy, and the specific heat are considered. We 
compare our results with experiments and investigate the 
infiuence of gravity. In Sec. IX we discuss dissipation and 
mutual friction for superfiuid helium below Tx- We show 
that ouji-approach reproduces the ansatz of Gorter and 
MellinkO for the mutual friction force and calculate the 
Gorter-Mellink coefficient A. The idea of the approach 
and part of the results, have been published already in a 
rapid communicationEj. 



II. THE MODEL AND FEYNMAN RULES 

Dynamic critical and heat transport phenomena in liq- 
uid ^He|-close to Tx are well described by model F which 
is givenlla by the Langevin equations for the order pa- 
rameter '0(r, t) and the entropy variable m(r, t): 
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where 



+ ^Xo^m^ +lom\'ip\^ - hom] (2.3) 

is the free energy functional and 00 and 9m are Gaussian 
stochastic forces which incorporate the fluctuations. The 
heat current Q is imposed by boundary conditions. The 
gravity is included via the temperature parameter to{z) 
in (2.3) which is related to Tx{z) and depends linearly 
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on the altitude z. Usually, the model is treated by field- 
theoretic means. The perturbation theory in terms of 
Feynman diagrams is genetated by the Janssen-De Do- 
minicis functional integralE3 



Z = 



J DipDijjDmDm cxp{J} (2.4) 



with the functional 
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Here '0 and to are auxiliary fielda^d which are needed for 
a proper construction of|4he perturbation theory. While 
in the previous theorieaj'cl only the free Green's func- 
tions were needed, here we will evaluate the self energy 
in leading order. For this reason, we briefly describe the 
Feynman rules by which the diagrams and the terms of 
the perturbation series are constructed. We decompose 
the fields into a mean-field and a fluctuating contribu- 
tion according to V' = V'mf + Sip and to = m^i + Sm 
where ipm i and TOmf are solutions of the model- equa- 
tions ( p^ll ) and (2^) without the stochastic forces. Here 
we assume a zero mean-field order parameter V-'mf — 
above and below T\, so that tp — H-'- The boundary 
conditions, which imply the heat current, require 



nf(z) = -q{xo/\Q)iz - Zo) 



(2.6) 



where q is the entropy current related to the heat current 
Q in physical units hy q = Q/ksTx. We decompose the 
functional J in powers of the fluctuating fields according 
to 



J = J2 + J3 + J 4 ■ 

The quadratic term is given by 

J2 = y d'^r J dt [Ao(Vto)' + r'o\ 



(2.7) 



with the differential operators 

Lo^dt+ ro[ro(z) - + 27oTO,„f(z)] 
-«(.9o/Xo)[nT-mf(2) - Xo^o] 

and 

No = dt- (Ao/xo)V2 . 
The third and fourth order terms are given by 



(2.8) 

(2.9) 
(2.10) 



J3 = j cf'r J dt{rn[XojQV\rr+9o'Vlm{rvr] 
-bsSmri^-blSmri)} (2.11) 

and 

J^^ [ d'^r [ dt{-birr^i^-Km^^ , (2.12) 



64 = 2roUo - igolo/"^ (2-13) 



respectively, where 
&3 = ro7o - jgo/2xo , 



are complex parameters. While the free Green's func- 
tions are obtained from the quadratic term J2, the in- 
teraction vertices are obtained from J3 and J4. In order 
to obtain a compact notation of the Feynman rules we 
combine the fields into vectors 



TO 

6m 



(2.14) 



where the index a = 1, 2 distinguishes between fields with 
and without tilde. Then the several contributions of the 
functional J can be written as 



J2 = - 



■J rft[^';ifo,a/3*/3 + ^MaEo,c.fiMp] , (2.15) 



J3 = 

J4 = -5 

Here 



d'^7 



d-^r 



dtB3.a0',Ma^*p^^ , 



2-^0 



(i?o,a/3) -y N+ 



(2.16) 
(2.17) 

(2.18) 
(2.19) 



are 2x2 matrices with differential operators as the el- 
ements. The free Green's functions are obtained by in- 
verting these matrices according to 

Go,a/5(r,i;r',t') = (^'„(r, t)«'^(r', <')>o 

^K-Xp5{v~v')5{t^t') , (2.20) 
7?o,a/3(r,i;r',t') = (M„(r, t)M/3(r', t'))o 

= ^^o:a/3'^(r-r')<5(t-t') . (2.21) 

The higher-rank tensors B^^afs-f and B4^af3,-fS, which de- 
scribe the interactions between the fields, are obtained by 
comparing (|]l|) and ( |2.17D with ( p.llD and ( |2.12| ), re- 
spectively. While i?3.a^-y contains differential operators, 
-64,0/3,75 is symmetrized with respect to interchange of 
the index pairs (a, P ) and (7, S). 

Now, from ( 2.15| )-( 2.21 ) the Feynman rules are ob- 
tained easily. In Fig. |1| the elements for constructing 
the Feynman diagrams are shown. The free ip-Held 
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Green's function Go,Qa(r, r', i') is identified by a di- 
rected solid line (Fig. |I|a). The free m-field Green's func- 
tion £)o,a/3(r, r', t') is identified by a dashed line (Fig. 
|f|b). The m-field interacts with the ij: fields by the three 
vertex shown in Fig. ^c. Furthermore, the fields inter- 
act with each other by the four vertex shown in Fig. 
For each three or four vertex an integration J cfir jdt 
must be performed. Further rules are applied as usual 
in field theory. The perturbation series is obtained as 
the sum of all possible Feynman diagrams which can be 
constructed from the elements shown in Fig. ^ 
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FIG. 1. The elements for constructing the Feynman di- 
agrams: (a) free 7/;-field Green's function, (b) free m-field 
Green's function, (c) three vertex, and (d) four vertex. 



III. THE UNRENORMALIZED GREEN'S 
FUNCTION IN HARTREE APPROXIMATION 

To obtain physical quantities we must first calculate 
the '0-field Green's function which is defined by 



G„^(r,i;r',t') = (^'a(r,t)vl'^(r',t')) 



(3.1) 



This Green's function can be expressed via the Dyson 
equation 



(3.2) 



in terms of the self energy I]Q^(r, i; r', t'). The pertur- 
bation series of the self energy is given by the sum of 
all irreducible Feynman diagrams with two amputated 
external solid lines which do not fall into pieces if any 
internal solid line is cut. A similar Dyson equation ex- 
ists also for the m-field Green's function Dai3{r,t;r' ,t'). 
However, for our calculations we do not need the latter 
Green's function explicitly. 

In order to account for effects beyond the perturbation 
theory we must resum the Feynman diagrams partially 
in an appropriate way. In the critical regime close to a 
second-order phase transition, infrared singularities oc- 
cur which must be resummed by renormalization and ap- 
plication of the renormalization-group theory. For model 



F the|B|G theory was elaborated up to two- loop order by 
Dohmt2l. All the renormalized coupling parameters, de- 
pending on a RG flow parameter were determinedlj by 
adjusting the superfluid density, the specific heat, and 
the thermal conductivity to the respective experimental 
data. Thus, model F can be used for explicit calcula- 
tions of physical quantities in the critical regime near T\ 
without any (further) adjustable parameters. 

Here we first apply an additional resummatiQii,"vstbich is 
well known in quantum many-particle physicsE3~Ej. We 
resum with respect to all self-energy subdiagrams so that 
the perturbation series becomes self consistent with re- 
spect to the ■(/'-field Green's function G. This means that 
now the solid lines are thick and identified by the exact 
Green's function G as shown in Fig. To avoid multiple 
counting of diagrams, only the irreducible diagrams are 
included in the perturbation series, which do not contain 
self-energy subdiagrams or equivalently which do not fall 
into pieces if any two of the internal thick solid lines are 
cut. fChe resummation was used first by Luttinger and 
WardEj and was formulated in terms ofpa-Legendre trans- 
formation by De Dominicis and Martiro. By truncating 
the self-consistent perturbation series of the self energjiX 
the conserving approximation of Baym and Kadanoffr^tr^l 
is obtained. 



(a) 
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= G„p(r,r;r', 







(b) L„B(r, f;r', r) = a_L p + „ 



r, t r,t 



P 

r, r r', r ' 



FIG. 2. (a) The exact i/i-field Green's function is identified 
by a thick solid directed line, (b) The self energy in Hartree 
approximation. 

Here, we approximate the self energy S by including 
only the tadpole diagrams as shown in Fig. ^b. This ap- 
proximation is equivalent to the Hartree approx ima tion 
in quantum mechanic£t3. The Dyson equation ( ^.2[ ) to- 
gether with the self energy in Fig. ^ are self-consistent 
equations which enable an explicit calculation of the un- 
renormalized Green's function G. However, since we con- 
sider liquid '*He in the critical regime near T\, a second 
resummation is necessary: the self-consistent perturba- 
tion series and hence the Hartree approximation must be 
modified by renormalization and application of the RG 
theory. 

We have several reasons to believe that the Hartree ap- 
proximation combined with the RG theory is successful 
for model F above and below T\ where the order parame- 
ter is always — 0. Moreover, we will show that the ap- 
proximation includes vortices and mutual friction so that 
it may be a possible approach to describe the dissipa±iQji 
in the superfluid state observed in the experimentsEnlj. 
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First of all, if we generalize model F by replacing the 
complex field ^ by a vector = {'ipi, . . . , V'n) of n com- 
plex fields, then it turns out that the Hartree approxi- 
mation is exact for the Green's function G in the limit 
n ^ oo. For each closed loop of thick solid lines there 
will be a factor n. If we rescale the coupling parameters 
according to uq ~ n~^, 70 ^ n~^/^, and go ^ n~^/^ so 
that B^^afSj ~ and B4_ai3^s ~ n~^, then only the 

tadpole diagrams shown in Fig. ^ will be nonzero in the 
limit n — > cxD. It is well known that models involving 
Ginzburg-Landau functionals like the t/)^ model can be 
solved exactly in this limit (see e.g. Ref. ^). The same 
is true also for model F. While in the limit n ^ 00 the 
RG theory is not needed because the Hartree approxi- 
mation is exact, in our case for n = 1 the RG theory is 
necessary to obtain the correct critical behavior of the 
physical quantities near Tx. 

We may expect that beyond the Hartree approxi- 
mation an 1/n expansion may yield the proper cor- 
rections. However, this is not true. By experience 
with quantupijjfield theory of many-particle systems with 
degeneraciesEjO we have found that proper corrections 
are given by the modified self-consistent random-phase 
approximation (modified SC-RPA), where the modifica- 
tion is a gauge transformation which implies bosoniza- 
tion for a proper treatment of the low energetic collec- 
tive excitations. The method has been invented for two- 
dimensional electron systeips in the regime of the frac- 
tional quantum Hall peffectcS and tested for simple ex- 
actly solvable modelscj. We have found that the modi- 
fied SC-RPA can describe the superfluid transition and 
Bose-Einstein condensation in interacting boson systems 
where the average order parameter is (ip) = [due to 
phase fluctuations above and below the transitioned. For 
this reason we believe that the modified SC-RPA com- 
bined with the RG theory will be successful in the present 
case for model F. 

The main feature of the modified SC-RPA is that it 
implies a nontrivial spectrum for the Green's function 
G while the Hartree approximation docs not. This fact 
is important in quantum many-particle physics because 
the spectra of the quasiparticles are important physical 
results. However, in the present case for critical phenom- 
ena and second-order phase transitions nontrivial spectra 
are not essential while the application of the RG theory 
is important. For this reason, the Hartree approximation 
combined with the RG theory should be sufficient for our 
purposes. Even though the modified SC-RPA combined 
with the RG theory would be desirable, we expect only 
some corrections while the calculations would be much 
more complicated. 



A. Evaluation of the Green's function 

Now, we evaluate the self energy E and determine the 
unrenormalized Green's function G in Hartree approxi- 



mation. The tadpole diagrams in Fig. Qb imply a self 
energy of the form 

S„^(r, t; r', t') = -AK^p S{r - r') S{t - t') (3.3) 

where AKap are the elements of a 2 x 2 matrix which 
depend on the space coordinate z but do not contain 
differential operators. Applying the Feynman rules of 
Sec. H we obtain the matrix 







where 

AL = 2b3No^[Xo-fo'V^ns + goVJ^] + 264ns . 
Here 

= ^G22(r,<;r,i) 
is a density related to the entropy and 

J, = (Im[7^*VV^]) = lim Im[VG22(r,t;r',i)] 



(3.4) 

(3.5) 
(3.6) 
(3.7) 



is the superfluid current dens ity. Nq^ is the invers e o f 
the differential operator (2.10). Since Tig and Jg in (3.5) 
depend only on z, we may drop the time derivative in 
N^^. Then, from (U) we obtain 



AL = 2(64 - Xo7o63)ns - 263(xogo/Ao)(V2)-i(VJs) 



(3.8) 



Inserting 63 and 64 of (2.13) we obtain 

h - Xoloh = 2ro(t2o - 5X070) = SFoMo , (3.9) 

where uq = uq — \xq1o is the effective coupling between 
the %l) fields in thermal equilibrium after the entropy field 
m has been integrated out (see Ref. |l^). Since the heat 
current Q flows in z direction, only the z component 
Js,z = Js of the superfluid current is nonzero so that 



Thus, from (|3.8[) we obtain 



AL = ToliuoHs - 2xo7o(5o/Ao) ^Js 
+«5o(5o/Ao) d^^Js ■ 



Now, we define 



L = Ln 



AL . 



Then from (2.£) and (3.11) we obtain 



L^dt + Toiniz) - V^] - z-^Aro(z) 



2xo7o 



(3.10) 

(3.11) 
(3.12) 

(3.13) 



where 
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ri{z) = To(z) + 2xo7o —mmiiz) - ^ 



(3.14) 



and 



Aro(z) = 2xo7o — m^f (z) - ho - ^ (3.15) 

are effective parameters. Furthermore we define the ma- 
trix 



and obtain 



0,a/3 



^ 

1 r + 







(3.16) 



(3.17) 



Now, the Green's function G in H artre e appr oxima tion 
is obtained easily. We find that (|3.12[) and ( [3.16 ) are 
equivalent to the Dyson equation ( |3.2[) . Thus, as a result 
we obtain 



G^p{v, t- r', t') = K-^ 5{v - r') 5{t - t') 



(3.18) 



where Kaf3 is given by ( ^.17 ) together with (3.13). 
Clearly, the Green's function G in Hartree approxima- 
tion has the same structure as the free Green's function 
Go where just the parameters have been replaced by ef- 
fective parameters. As a consequence, the following cal- 
culations are considerably simplified because we may re- 
strict the considerations to the effective parameters ri{z) 
and /S.ro{z) only. We will derive self-consistent equations 
for ri{z) and Aro(z) to determine the effective parame- 
ters. Eventually, the Green's function G is obtained from 
( |3ll) together with ( ^Isj ) and 



B. Evaluation of and Js 

Next w e eva luate rtg a nd Js by inserting the Green's 
function ( 3. IS ) into ( ) and ( 3/7 ) . In Ref . ^ and Jg 



were evaluated for the free Green's function so that here 
we need not perform these calculations once again. Since 
the Green's function in Hartree approximation has the 
same structure, we may use the results of Ref. || with a 
slight modification. We just need to replace the free pa- 
rameters by the effective parameters ri{z) and Aro(2:) 
appropriately^. The basic assumption of the previous 
calculationgj was that the parameters ri{z) and Aro(z) 
in the Green's function are linear functions of z. For 



t he fr ee Green's function, where rig and Js in (3.14) and 
( 3.15 ) are omi tted, this is indeed true, because mmf(^) 
defined in (2.6) and to(z), which is related to T\{z), are 
linear functions of z. Thus, in the present case we must 
assume as an approximation, that ri{z) and Aro(2;) are 
linearized locally so that only the slopes r'^ and Arp are 
included but the curvatures and the higher-ord er de riva- 
tives are neglected. We will later show in Sec. VH that 
this assumption is justified. 



Now, from (3.22) and (3.25) in the second paper of 
Ref. ^ wc obtain 



n., = 2 $_ 



l+e 



,2{X) 



d'^p 1 
{2ttY ri(z) 



(3.19) 



and 
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Ar'o 



2r[, 2x070 



d'^p 1 
(2^ [ri(z)+p2]2 



(3.20) 



Here e is related to the space dimension d by e = 4 — d. 
The function ^^(X) is defined by the asymptotic series 



^ r(a + 3N) X 



T{a) 



N 



-TTT (3-21) 



and contains all the effects beyond linear response the- 
ory. In the integrals the parameter fo{z) of Ref. || has 
been replaced by the effective parameter ri{z). The di- 
mensionless parameter X is given here by 



X 



12[ri(z)]3 



„'2 



r" 



nv4xo7on 



where 



r'l = dzri{z) 



' 90 ^, 
^4xo7on ° 



Ar'o = 9.Aro(z) 



(3.22) 



(3.23) 



The entropy current q in the formulas of Ref. |^ is here 
replaced by or Arp times — Ao/2xo7o- We note that 
our formulas reduce to those of Ref. if we omit 71= and 



Js in the effective parameters ( ^.14P "^and (3T£) ^pd if 
we neglect gravity. The integrals can be evaluatedclilj in 
dimensional regularization so that we obtain 



ns{z) 



<5-i+./2(^) [ri{z)Y-''^ , (3.24) 



and 



1 I \ So 



Ar 



2r(, 2x070 e 



1-1 
2 



^,/2ix)[n{z)] 



-e/2 



(3.25) 



where = ^<jr(l - e/2)r(l + e/2), - f^d/(27r)^ and 
Qd = 2n'^^^ /T{d/2) is the surface of the d dimensional 
unit sphere. Clearly, ns{z) and Js{z) depend on z im- 
plicitly via ri{z) and the derivatives r[ and A^q. 



C. Self-consistent equations for 
the effective parameters 



Eqs. ( |3.14| ) and ( |3.15| ) together with ( |3.24| ) and ( |3.25D 

are self-consistent equations for the effective parameters 
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ri(z) and ArQ^z). The structure of these equations can 
be simphfied. First of aU we note that the effective pa- 
rameters and its derivatives can be related to the tem- 
peratures T{z), T\{z), and the heat current Q. We find 



Aro(z) 



/SH\ T{z) - To 
2xo7o( J^/ = 2xo7o 



(3.26) 



To prove this relation {SH/6m) can be evaluated explic - 
itly in Hartree approximation and comp ared with (3.f 5). 
Equivalently, we take the average of (2^) and obtain 
dt{m) -f Vq = where q 



q = -AoV 



\ dm I 



.9oJs 



(3.27) 



is the entropy current. Since z is the only space coordi- 
nate, this equation can be rewritten in the form 



d /5H\ 
dz \ 5m I 



d_ 
9i 



— mmf (z) 
Xo 



(3.28) 



Thus , integrating this equation and c ompa ring with 
( 3.15 ) we obtain ( 3.26| ). We note that ( 3.15| ) contains 
integration constants via m-^^iz) and d~^Js- Since m is 
the entropy density divided by fee, the quantity (SH/dm) 
is a temperature difference divi ded b y Tx- This fact ex- 
plains the last equality sign in ( 3.26 ). Here Tq is a con- 
stant reference temperature which may be arbitrary. In 
the denominator the z dependence of T\ due to gravity 
is very small and may be neglected. 

The temperature parameter ro{z) is defined by 



roiz) = To{z) + 2xo7o(/io + (^)) 
= To{z) + 2xo7o/io + Ara{z) . 



(3.29) 



In thermal equilibrium this parameter is the coefficient 
of in the energy functionalr-after the entropy vari- 
able m has been integrated outEj. It is related to the 
temperature by 



roiz) - roc = 2xo7o 



T{z)-Tx{z) 



= 2xo7o 



Ar(z) 

Tx 



(3.30) 



The critical value is tqc = in one-loop approximationll3 
and hence also in Hartree approximation. We find that 
the first line on the r ight- hand side of ( 3.14| ) is identified 
by ro{z). Thus, Eq. (3.14) simplifies into 



ri{z) = ro{z) + AuqUs 



(3.31) 



We resolve this equation with respect to ^0(2:) and insert 
(3.24) for rig. Then we obtain 



roiz) = niz){l + ^ A,a>_i+,/2(X) [ri(z)]-/2} 



(3.32) 



Next, we take the derivative of this equation with respect 
to z and obtain 



r'o = r[{ 



8uq 



iX)[r,iz)r^^'} 



Furthermore, the derivative of ( 3.15 ) yields 



= -2xo7o 



- + -Js 
Ao Ao 



(3.33) 



(3.34) 



Resol ving this equation with respect to q and inserting 
( 3.25 ) for Js we obtain 



Q 



kuTx 



-Ao 1 



.9§ 1 



2Aor' e 



-A, 1-^ 



c$,/2(X)[ri(z)]-/2}^ . (3.35) 



On the other hand from (3.26) we find that Ar^ is related 
to the temperature gradient dzT by 



Ar'o = 2xoloTx^d,T 



Finally, Eq. (|3.29D implies 

r'^ = Ar[, - 2xoloTx^d,Tx , 



(3.36) 



(3.37) 



where we have identified Tq = ~2xa^QT^'^dzTx- Clearly, 
the difference between Tq and Arg is due to the gradi- 
ent of Txiz) which is the effect of gravity. Thus, in a 
microgravity environment Tq and Ar^ are equal. 

Now, the self-consistent equations which allow the de- 
termination of all the effective parameters are given by 
(|332|), (^), and (|3|) together with ( |3^ ), ( pO| ), 
( 3.37 ), and ( 3.22 ). These are seven equations for seven 
unknown variables ^0(2:), Aro(2;), ri(z), Tq, Arp, r'l, and 
X. As an input we need the temperatures T(z), Txiz)^ 
the difference Ar(z) = r(z) - Ta(z), the heat current 
Q — k^Txq, and the gradient d^Tx for a given space 
variable z. Since the seven equations do not depend ex- 
plicitly on 2, we do not need the temperature profiles as 
functions of z. Instead, we obtain the temperature gradi- 
ent dzT from (3.36) so that the temperature profile r(z) 
can be calculated by integration. Eventually, for given 
AT and Q physical quantities like the specific heat and 
the thermal conductivity can be calculated. This will be 
done in Sees. ^ and VIII. 



IV. RENORMALIZATION AND APPLICATION 
OF THE RENORMALIZATION-GROUP THEORY 



In ( 3.32| ), (3.33), and ( [3.35| ) the first-order terms ex- 
hibit infrared divergences at criticality where ri(z) — > 
while the function <I>q(X) is of order unity. For this rea- 
son the renormalization of these equations and the ap- 
plication of the RG theory are necessary to achieve a 
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resumniation of the infrared divergences and a proper 
treatment of the critical fluctuations. We use the con- 
cept of renormahzation by minimal subtraction of dimen- 
sional poles. The calculations are performed at fixed di- 
mension d = 4 — e (i.e. no e expansion is applied). For 
model F this renormahzation scheme is described in Ref. 
[l^ . The renormahzation factors of the fields V' and -0 
are = = 1 in one-loop order and hence also in 
Hartree approximation. Thus, the Green's function G is 
not renormalized here. As a consequence, also the oper- 
ator L defined in ( [3.13 ) is not renornaalized. The param- 
eter Fq is renormalized according toE^ Fq = Zf^T where, 
however, in one-loop order and in Hartree approximation 
Zr = 1. Thus, from ( [3.13| ) we conclude that ri{z) is not 
renormalized. 

j^he parameters 70 and go are renormalized according 



t. 



X0l0=j{X0Zrn)'/'Zr{fiyAa)'/^ 



(4.1) 
(4.2) 



Furthermore, we renormalize Aro(z) = ZrAr(z). Conse- 
quently, we find 



30 



2xo7o 



Aro(z) 



_9_ 

27 



Ar{z) 



(4.3) 



which is cons istent with the requirement that the last 
term in (3.13) is not renormalized. Thus, in terms of the 
renormalized parameters the operator ( 3.13| ) reads 



L^dt+ F[ri(z) - V2] - i{g/2j)Ar{z) 



(4.4) 



Next we renjOsmalize ( 3.32 ). For this purpose we need 
the relationsEJ 



uo = uZuZ^ (^V^d) 



(4.5) 
(4.6) 



where tqc = and Z^ = 1 in Hartree approximation. 
We sep arate the ultraviolet divergence on the right-hand 
side of (3.32) which here in dimensional regularization is 
a pole ~ 1/e. By choosing 



Zr = Zu = 1/[1 - 8u/e] 



(4.7) 



the ultraviolet divergence is canceled. Eventually we ob- 
tain 



"(z) = ri(z)|] 



8u 
e 



l+e/2l 



(x,(iiM)-"_,]} 



/2 



(4.8) 



Analogously we separate the ultr aviol et div ergence on the 
right-hand side of ( |3.33| ). Using (m)-(|^ we obtain 



ri(z)Nj-<^/2 

7' 



(4.9) 



We replace t he pa rameters in the dimensionless variable 
X defined in ( 3.22| ) by the renormalized ones. It turns out 
that all Z factors cancel so that X is not renormalized. 
Thus, we obtain 



X 



1 



12[ri(z)]3 



r \4:jr J 1 V47r' 



(4.10) 



For convenience we replace the reaormalized couplings 
by the dimensionless combinational^ w = F/A, F — g/X, 
and / = F'^/w'. Furthermore, we introduce the dimen- 
sionless effective parameters 



p = 




P' 




(4.11) 




Ar(z)/M2 , 


Ap' 


= Ar'/n^ , 


(4.12) 


Pi = 


ri(z)/^2 , 


P'l 




(4.13) 



Since the coordinate z does not appear explicitly in the 
self-consistent equations for the effective parameter s, w e 
omit z as an argument from now on. Then, Eqs. ( |4.8| )- 
(4.10) can be rewritten as 



P = 

X ' 



8u 



e 

8u 



12p^ 



P? 



/,{X)pr^'-l\} , 

i-|) <p,/,ix)pr^'-i\} 

+ 2^" ( ^ Ap')p' ( ^ 



respectively. ^ 
The entropy current is renormalized hya 



(4.14) 

(4.15) 
2n 



(4.16) 



(4.17) 



This equation together with (4.2) implies that the ratio 

Q/gokBTx - q/go = (g'-^V.g) (Ad/p')'/' (4.18) 
need not be renormalized because the Z factors cancel. 



In (3.35) we separate the ultraviolet divergence and re 
place the coupling parameters by the renormalized pa 
rameters. Additionally, we need the renormalizationt 
Aq = xqZ^^X- Using the Z factor product 



Z„,Za = 1/[1 - f/2e] 



(4.19) 



we find that the ultraviolet divergence is canceled. Then, 
in terms of the renormalized couplings and effective pa- 
rameters we obtain 



'P 



£-3 



Ad 
27F 

1 - 



{'4 



<i>,/2(X)pr/'-l]}Ap' (4.20) 



where the left-hand side need not be ren orma liz ed be - 
cause of ( 4.18|). Finally, we renormalize (3.26), (3.30) 
(3.36), and ( |3.37 ) and obtain 
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Ap = t-1[T(z)-To]/Ta , 

p = T-\T{z)-Tx{z)]lTx , 
Ap' =T-HfiTx)-^d,T , 
p' = Ap' -T-HiiTx)-^d,Tx 



respectively, where 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



(4.25) 



In these equations Zm does not cancel. We note that 
the renormalization is exact in all the above equations. 
This means that we need not expand the Z factors in 
powers of the renormalized couplings. The reason of this 
fact is that the Hartree approximation is exact in the 
limit n — > 00 for model F with an n-component com- 
plex order parameter. On the other hand, the 2_factors 
do not agree with those of the previous theoriealj in one 
loop-order because the Hartree approximation is not a 
loop expansion. (The correct one-loop Z factors would 
be obtained if we would consider a Hartree-Fock approx- 
imation and include both the Hartree and the Fock term 
in the self energy S.) 

By the renormalization a characteristic length scale is 
introduced which is described by the parameter p,. The 
RG theory is based on the fact that this length scale 
is arbitrary and may be changed according to p pi, 
where I is the RG flow parameter. As a consequence, the 
renormalized couphng parameters u{l), j{l), w{l), F{1), 
and f{l) depend on I. Furthermore, also the Z factors 
depend on I. Now, the dimensionless parameter defined 
in (4.25) reads 



\\[\Zrr,{l)' 



/2 1 



Xo^m(0^ 27(/) 



(4.26) 



For convenience we will use r as the RG flow parameter 
instead of I because r is closely related to the reduced 
temperature by ( [1.22 ) and the renormalized coupling pa- 
rameters It [r] , 7 [r] , w[t\, F\t\ , and / [t] were determined 
as functions of r in Ref. O. We identify pi = by 
the correlation length ^ = ^(r), which in the asymp- 
totic region is given by ^(r) = ^o'''"''- The identification 
pi = is correet in one- loop order, corrections appear 
in higher ordersEJ. 

Now, we write the self-consistent equations for the ef- 
fective parameters in a form which is appropriate for the 
numerical evaluation. For this purpose we eliminate some 
of the dimensionless parameters and introduce some new 
parameters. First of all we note that the asymptotic se- 
ries ( p.2_l[ ) is not useful to evaluate the function <I>q(X). 
In Rcf . an integral representation was found by 



where C = (—X) and 



•^a(C) = / dww""^exp(- 

"'0 



vO . 



(4.27) 



(4.28) 



The integral is well defined for a > 0. For a < the 
function J-'a{C) is obtained by analytical con tinua tion in 
a or equivalently by partial integration in (4.28) to re- 
move the ultraviolet divergence. Now, we introduce the 
parameter 



1 

12 



■■[T\ / F[t]Ap' \ _ / ^irj Ap- \ 



F[t]Ap' \2 



w'[t] V 47 [r]w'[r 



47 [t]w' [t 

(4.29) 



so that X = —a/p\ or equivalently pi = a^^^C,. In the 
following we eliminate pi and X in favor of C and a. For 
convenience we define the amplitudes 



A = e-i[$_i+,/2(^)pr -1] 
_ ir o-^Ai 

~7lr(-i + e/2)'' -^-1+^/2 

A,=e-'[{l^e/2)<^>,/,{X)pr^' 



(C) - 1 
-1] 



r(-l + e/2) 



Then we rewrite (4.14) as 

p = cri/^C{l + 8u[r]^} 



(4.30) 
(4.31) 
(4.32) 
(4.33) 

(4.34) 



To eliminate p we insert this into ( [4.22| ) . Resolving with 
respect to the temperature difference we obtain 



AT{z) = T{z) - Tx{z) = Txra^/^C {l + 8u[t]A} 



(4.35) 



Next we resolve (4.20) with respect to Ap' and obtain 

''"'^ '{l + (/[r]/2)Ai} . 

(4.36) 



^ _27[t]F[t] / 



Resolving (4.15) with respect to pi and eliminating p' by 
inserting (4.24) we obtain 



Pi 



Ap' 



l^dTx 
tTx dz 



{1 + 8u[t]Ai} 



(4.37) 



Finally, from ( 4.23| ) we obtain the temperature gradient 

d,T = Tx {t/0 Ap' . (4.38) 

Until now the RG flow parameter r is arbitrary. We 
must choose r so that an optimum resummation of the in- 
frared divergences in the perturbation series is achieved. 
From our experience we find that the condition 



CT^/^ {8 + C-16u[t]AC) ^ 1 



(4.39) 



is an optimum choice for fixing the RG fl ow p arameter t. 
For Q — > in thermal equilibrium Eq. (4.3!;) reduces to 
the well known flow parameter conditions of the previous 
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theoriesEi above and below T\ . The integral ( 4.28 ) can be 
evaluated asymptotically for large positive and negative 
Q. For C ^ +1 we find 



^a(C)«r(a)r" 



(4.40) 



which implies <i>Q,(X) « 1. Consequently, from ( 4.30| )- 
(4.33) we obtain the amplitudes 



A 
A, 



e/2 



-M(i-6/2)pr^^-i] 



(4.41) 
(4.42) 



The flow parameter equation ( 4.39| ) reduces to pi(l — 
16u[r]^) = 1 which implies 



r(0/(/iO'-P-Pi = l 



(4.43) 



Eq. ( |4.22| ) imphes T > Ta and r = (T - Tx)lTx so that 
the RG flow parameter r is identified by the reduced 
Indeed, Eq. (|4.4^ 



temperature. 

equation of the previous theoriesi 
for T > Tx- On the other hand for ( < 



is the flow-parameter 
in thermal equilibrium 
1 we find 



T^iC) « (7r/3)i/2(_c/3)"/2-3/4exp{2(-C/3)'/n (4-44) 

which is exponentially large. Consequently, the ampli- 
tudes A and Ai are exponentially lar ge so that in ( [4.34 ) 
and in the flow parameter condition ( 4.39 ) only the last 
terms are relevant. Eliminating A we obtain 



2r{l)/{pilf = -2p= 1 



(4.45) 



which is the flow-para meter equation of the previous the- 
ories for T <Tx. Eq. ( |I2^ ) implies t ^ ~2{T - Tx)/Tx 
and T < Tx- We conclude that our preseajt theory for 
nonzero Q reduces to the previous theoriest2l for C S> +1 
in the normal fluid region well above Tx and for C ^ ~ 1 
in the superfluid region well below Tx- 

We have derived seven equati ons given by ( 4.29 ), 



( |OlD , dog ), (p5[)-(p7|), and ( |4.39| ) which we have 
published alrcadp^ for d — i and e = 1 in a rapid 
communicatiorO. These equations contain seven vari- 
ables C: cr, T, A, Ai, Ap' , p'l which can be determined 
uniquely by solving the equations supposed the temper- 
ature difference AT = T — Tx and the heat current Q are 
known. The remaining effective parameters, which we 
have eliminated, can be determined afterwards. In prac- 
tice we have solved the equations for d — 3 dimensions 
and e = 1 in the following way. While the heat current Q 
is assumed to be constant we take as a variable which 
we vary in the whol e int erval —oo<^< +oo to scan all 
temperatures. Eq. ( 4.39 ) is solved explicitly to obtain a 
as a function of C and r. Then we solve the equations 
numerically by adjusting the flow parameter r and even- 
tually determine the temperature difference A T and the 
temperature gradient dzT as functions of ( by ( 4.35 ) and 
(p8|). 



As an input we need the dimensionless renormalized 
couplings m[t], 7[r], w[t] = w'[t] +iw"[T], F[t], and /[r] 



as functions of r which have been determined by DohmEl. 
Furthermore, we need the parameter go which is related 
to the entropy at Tx- For liquid, helium at saturated va- 
por pressure this parameter i£d go = 2.164 x 10^^ s~^. 
To calculate the correlation length ^(r) = ^o'''"'^ as a 
function of r we use the exponent v = 0.671 and the am- 
plitude ^0 = 1.45 X 10~* cm which were determined ex- 
perimentally in Refs. 31 and |3^. There arc no adjustable 
parameters. 



V. THERMAL CONDUCTIVITY 

A. Numerical evaluation of the thermal conductivity 
and comparison with experiments 



We eliminate Ap' from ( 4.36| ) and ( 4.38| ) and resolve 
the resulting equation with respect to Q. Then we obtain 
the heat transport equation 



Q 



-At^T 



where 



At — 



gokBAd {1 + (/[t]/2)Ai} 



27[t]T[t] 



(5.1) 



(5.2) 



is the thermal conductivity. Inserting the dimensionless 
parameters into (5.2), which we calculate for given AT 
and Q as described above by solving the seven equations, 
we obtain the thermal conductivity At = At(AT, Q) as 
a function of AT and Q. The result is obtained without 
adjustable parameters. We plot the thermal resistivity 
Pt = 1/ At logarithmically as a function of AT = T — Tx 
for given heat currents Q. In Fig. ^ our result is shown 
for Q = 42.9 /xW/cm^ as solid line. First of all we flnd 
that Pt is nonzero and At is finite for all temperatures 
above and below Tx- For T well above Tx and well below 
Tx we approximately find asymptotic power laws, which 
we will discuss in the next subsection. 

In the intermediate region, where T is close to Tx, the 
curve interpolates the two approximate power laws and 
shows a point of maximum slope. We may interpret this 
point as the superfluid transition and deflne the related 
shift of the critical temperature ATx{Q), which irt-Fig. 
1^ is indicated by the arrow. In the previous theorytil the 
formula 



ATa(Q) = -MT; 



V^ofeTA/ 



(5.3) 



was derived for the shift of the critical temperature with 
the exponent x ~ [{d — 1)H~^ = 0.745 and the constant 
M = 2.90 for d = 3 dimensions. Here we use ( |5.3|) as a 
flt formula for the point of maximum slope in Fig. H. By 
varying the heat current Q we flnd nearly the same expo- 
nent X = 0.745 where deviations due to nonasymptotic 
effects of the dynamic RG theory are very small here. 
Furthermore, we flnd the constant M = 3.17 which also 
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is nearly the same. Thus, the point of maximum slope in 
Fig. ^ may indeed be identified as the superfluid -transi- 
tion. However, in contrast to the previous theoryo, here 
AT\{Q) and the superfluid transition are not sharply de- 
fined for nonzero Q because the curves of the physical 
quantities are smooth. 




-8 -4 

AT=T-T^ [^K] 

FIG. 3. The thermal resistivity pT = 1/At logarithmi- 
cally as a function of AT = T — T\ for the heat current 
Q = 42.9 piW/cm'^. Tl*e solid line represents our theory. The 
data of Liu awd Ahlersa are shown as crosses, and the data of 
Baddar et alJiJ (fit formula) are shown as dashed line. The 
arrow indicates l\T\{Q). The dash-dotted line extrapolating 
the crosses indicates ATc((3). For comparison, the thermal 
resistivity for (5 = (theory) is shown as dotted line. 

For zero heat current (linear response limit) the ther- 
mal conductivity is shown as dotted line in Fig. |^, which 
is finite for AT > 0, diverges at AT = 0, and is infinite 
for AT < 0. Clearly, the solid line of our present theory 
(At (AT, Q) at nonzero Q) approaches the dotted line (At 
at Q = 0) asymptotically in the normal fluid region for 
large positive AT. Furthermore, the present theory re- 
produces the previous theory of Ref. ^ for AT > 0. Some 
small deviations of the two approaches from each other 
occur because here and in Ref. || the RG flow parame- 
terr is determined by different conditions (the condition 
.39 ) here differs from (4.49) in the second paper of Ref. 



M). However, these deviations are within the errors of the 
RG theory and heiice are not serious. 

Liu and AhlersEl considered a vertical heat flow up- 
wards in liquid "^He. They measured the temperatures 
at the bottom and at the top of the cell and deter- 
mined the thermal conductivity At(AT, Q) in two dif- 
ferent ways. First, they used a fit formula with power 
laws for At (AT, Q) and adjusted the exponents and am- 
plitudes. Secondly, they obtained At by the differential 
formula (24) of the first paper of Ref. I In Fig. I the data 
obtained by the differential formula are shown as crosses. 
For lower temperatures the data are extrapolated by the 
fit formula shown as dash-dotted fine. While AT\{Q) is 



indicated by the arrow, the data show the second tran- 
sition temperature ATdQ) where the dash-dotted line 
drops down nearly vertically. In the normal fluid region 
for AT > ATx{Q) the experimental data agree with our 
theoretical prediction (solid line) within the accuracies 
of theory and experiment. However, in the superfluid. re- 
gion for AT < ATx{Q) the data of Liu and AhlersB do 
not agree with our theory. For temperatures AT in the 
interval AX'c(Q) < AT < ATa(Q), the so called dissipa- 
tive regioriQ, the experiment finds a much larger thermal 
resistivity pT = 1/-^t, by about a factor of 20, than our 
theory predicts. Then suddenly, when the temperature 
AT approaches and drops below the second transition 
temperature ATc((3), the experimentally observed ther- 
mal resistivity is so sraall, that it was not detected any 
more in the experimentu. Our theory does not predict the 
second transition tempetature ATc((5), which has been 
found in the experimentsB'ErEl. While the solid line in Fig. 
1^ has a inflection point with maximum slope at AT\{Q), 
nothing unusual is found at ATc(Q). 

The experiments of Refs. ^ and ^ were performed by 
measuring the temperatures at the bottom plate and at 
the top plate, where the heat current flows into and out 
off the helium. For this reason, these experiments may be 
influenced considerably by surface effects. In superfluid 
"^He the dissipation of the heat current is caused by cre- 
ation of vortices, where the thermal resistivity pT = 1/At 
is proportional to the density of vortices in the helium. 
Near the bottom and top surfaces additional vortices may 
be created which enhance the vortex density there. This 
effect may possibly explain|-|the strong enhancement of 
the experimentally observedD thermal resistivity for tem- 
peratures AT in the interval ATc(Q) < AT < ATx{Q) 
(see crosses in Fig. ^). ■— ■ 

In a recent experiment Baddar et al.ll3 measured the 
temperature gradient dzT in superfluid helium for sev- 
eral heat currents Q. To exclude surface effects at 
the bottom and top plates the temperatures were mea- 
sured by sidewall thermometers only. The thermal con- 
ductivity At(AT, Q) wa s then obtained from the heat 
transport equation (5.1). A power-law fit formula for 
At (AT, Q) was found which is valid for a wide range of 
heat currents Q and temperatures AT sufficiently below 
ATa(Q). In Fig. I the thermal resistivity pT = 1/At 
for Q — 42.9 /xW/cm^ represented by this fit formula is 
shown as dashed line. Clearly, the temperature depen- 
dence of the experimental data agrees qualitatively with 
the theoretical prediction (solid line). However, the abso- 
lute values of the measured thermal resistivity are about 
a factor of 20 smaller than the theoretically predicted 
values. While the new experiment of Baddar et alO is 
believed to be a better and more direct measurement of 
the thermal conductivity or resistivity, there remains a 
disagreement between experiment and theory. 

We may possibly explain the discrepancy in the fol- 
lowing way. In superfluid "^Hc a honiogeneous heat cur- 
rent Q represents a metastable statecl. For the creation 
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of vortices energy barriers must be overcome. This fact 
keeps the rate of vortex creation low so that the vortex 
density and hence the thermal resistivity are small. On 
the other hand our theory is based on the approxima- 
tion where the complex order parameter is replaced by 
a vector ^' = {ipi, . . . ^tpn) of n complex components in 
the limit n — s- oo. In this limit the heat current Q is al- 
ways unstable so that the rate of vortex creation is higher. 
Consequently, in our theory the vortex density and hence 
the thermal resistivity are expected to be larger. How- 
ever, the large discrepancies in Fig. ^ indicate that the 
vortex density is a very sensitive quantity which may be 
influenced strongly by the kind of the approximation in 
theory and by certain conditions in the experiment. 



B. Asymptotic formulas for the thermal conductivity 



For temperatures T well above Tx Eqs. ( [4.4lD -( ^.43D 
imply A = Q and Ai = —1/2 so that At depends only 
on T = AT/T\ but not on the heat current. Thus, well 
above Tx the heat transport described by ( |5.lD is linear. 
From (^^) we recover the well known result of Ref. |l3| for 
At at infinitesimal Q in one-loop order. Asymptotically 
in leading order we find 



-(d-2) 



AT 



(5.4) 



which diverges in the limit AT — > 0. For cZ = 3 includ- 
ing nonasymptotic effects of the dynamic RG theory one 
obtains At ~ /\j^-OAi^ Thus, near Tx the thermal con- 
ductivity is strongly enhanced by critical fluctuations. 

However, at fixed nonzero Q for temperatures T close 
to Tx the flow parameter r will reach a minimum value 
and the correlation length ^ will reach a maximum value, 
so that At remains finite even at Tx ■ This fact was found 
previously in Ref. |^ and is seen clearly in Fig. ^. The ther- 
mal conductivity At = At (AT, Q) becomes Q dependent 
so that the heat transport is nonlinear. The crossover 
from linear to nonlinear heat transport happens for tem- 
peratures AT belowd 



AT„i(Q) =M„iTAf 



ygoksTx 



(5.5) 



with the exponent x ^ [{d — l)i^]^^, where M^i is a con- 
stant of order unity. Eor d — 3 the values x — 0.745 and 
Mni w 2.8 were foundB. _ 

While the previous theoryS is valid only for AT > 0, 
the present theory works also for lower temperatures in 
the superfluid region. Well below Tx the thermal conduc- 
tivity (5.2) can be evaluated asymptot ically . For C, ^ —5 
the function J-a{C) approximated by (4.44) is exponen- 
tially large. Consequently, t he am plitud es A and Ai are 
exponentially large. From ( 4.31 ) and ( 4.33 ) we obtain 
the ratio 



Ai 
A 



>_l+,/2(C) 



(5.6) 



Eqs. ( |436| ) and ( |437| ) reduce to 



F[t] Ap' 



47(7 



Ad V^o^bTa/ Ai 
1 ^ dTx I 1 
T Tx dz 8u[t 



(5.7) 
(5.8) 



We assume that the heat current Q is suffici entl y large so 



that gravity effects may be neglected and (5.8) is much 



smaller than (5.7). Then, from (4.29) we obtain 



\1A\ \gokBTx) A\ 



(5.9) 



On the other hand, the flow parameter equation ( 4.39 ) 
reduces to 

cri/3i6u[r]^(-C) = 1 . (5.10) 

Now, Eqs. (5_^), ( |5.9| ), and ( ^.10| ) are three equations for 
A, Ai, and a. Eliminating a we obtain the amplitudes 



Ak- 



^d 



■A i8M[r])3 I, Q^d-1 



2^/3 (8«[r])3VQe 



(5.11) 
(5.12) 



Eventually, from (| 
ity 



3) we obtain the thermal conductiv- 



At 



go fee F\t\ 



(5.13) 



where r = — 2AT/Tx. The variable ( depends only 
weakly on r, i.e. logarithmically. Thus, asymptotically 
in leading order we find 



At ~ T-^^'^-^'^Q-^ - (-AT)(3''-'*)''-iQ-2 (5 ^^4) 

For d = 3 including nonasymptotic effects we obtain 
At (-AT)2-4Q-2. 

In their experiment Baddar et al.EII have found that for 
a wide range of heat currents Q the thermal conductivity 
can be expressed in terms of the power-law formula 



At,cxp = Ao[(-AT/Ta)(Q/Qo)-°-'"' 

~ (-AT)2 «Q-2-53 



2.8 



(5.15) 
(5.16) 



where Ao = 1 W/cmK and Qa = 393 W/cm^ Clearly, 
this power-law form ula h as the same structure as the 
asymptotic formula ( 5.14 ) of our theoretical prediction. 
However, the two formulas do not agree quantitatively 
with each other. The exponents of the power laws differ 
by about 20%. Furthermore, the amplitudes differ by a 
factor of 20 which means that At,oxp is about 20 times 
larger than the theoretical At (see dashed and solid line 
in Fig. ^). The discrepancies may possibly be due to the 
approximation we have used in our theory. Further theo- 
retical and experimental work is necessary to clarify the 
origin of the discrepancies. 
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C. Influence of gravity 

On earth gravity imphes a spatiaUy dependent super- 
fluid transition temperature T\{z) with a nonzero gradi- 
ent dzT\ — ±1.273 ^K/cm. Supposed the heat current Q 
is flowing in the z direction the gradient dzT\ is positive 
for heat current upwards and negative for heat current 
downwards. In our theory the gradient d^Tx is incorpo- 
rated in (4.24) or equivalently in (4.37). We find that on 
earth the effects of heat current Q and of gravity have 
equal magnitude for Q « 65 nW/cm^. For larger Q the 
heat current is dominating while for smaller Q gravity is 
dominating. We find that for Q > 1 /iW/cm^ the effects 
of gravity are very small so that in this case gravity can 
be neglected. The experiments by Liu and AhlersO, Mur- 
phAL,and MeyerQ, and the new experiment by Baddar et 
al.EZl are performed at heat currents Q which satisfy this 
condition. 

We have performed our calculations for a positive and 
a negative gradient dzT\ representing gravity on earth 
and also for a zero gradient dzT\ which corresponds to 
a microgravity environment in space. As a result, for 
heat currents Q 100 nW/cm^ the thermal conductivity 
At (AT, Q) is nearly independent of gravity. On the other 
hand, for heat currents Q ~ 80 nW/cm o ur th eory fails 
if gravity is present, because a defined in ( 4.29 ) changes 
sign so that no solution of the seven equations for the 
dimensionless-«arameters in Sec. IV can be found. 

Day et al.E3 investigated the superfluid-normal-fluid 
interface in "^He and measured the thermal conductiv- 
ity At for very small heat currents Q in the inter- 
val 20 nW/cm < Q < 6 /xW/cm^. The heat cur- 
rent flows upwards so that the gradient dzT\ is pos- 
itive. Clearly, this experiment explores the crossover 
from the gravity dominated region to the heat current 
dominated region. For heat currents Q > 100 nW/cm^ 
the experimental data agree quite well with our theo- 
retical prediction for At. This has Jbeen demonstrated 
in our previous rapid communicationEfl for temperatures 
AT > ATxiQ) (see Fig. 2 therein). In the superfluid 
region for AT < ATx{Q) the thermal resistivity is very 
small so that here experimental data are not available 
with sufficient accuracy on a logarithmic scale. 

For very low h&at currents Q < 40 nW/cm^ the 
data of Day et al.E3 indicate that the thermal resistiv- 
ity Pt — 1/At(AT, (5) is a smooth function of AT in 
the limit Q if gravity is present. This fact means 
that gravity prevents the system from reaching the criti- 
cal point of the superfluid transition, so that all physical 
quantities are smooth and nonsingular near T\ . Unfortu- 
nately, our theory fails for these low heat currents in the 
gravity dominated region. However, an alternative at 
proach is possible which is equivalent to Onuki's theorj 



In m ean- field approximation the model-F equations ( 2.1 ) 
and ( |2.2| ) can be solved numerically as partial differen- 
tial equations to obtain the order-parameter profile and 
the temperature profile of the superfiuid-normal-fluid in- 



terface. To include the critical fluctuations the model- 
F equations are renormalizcd and the RG theory is ap- 
plied. The approach is a renormalizcd mean-field theory 
which works for all heat currents Q and gravity, even for 
very small Q in the gravity dominated region. In this 
way we obtain a thermal resistivity pT = 1/At(AT, Q) 
which agrees qualitatively with the experimental obser- 
vation of Day et al.E3: it is a smooth function of AT and 
nearly independent of Q for Q ^ 10 nW/cm^ including 
the limit Q ^ 0. In agreement with the experimentE^ we 
find that STg « 15 nK is the relevant temperature scale 
on which the critical singularity is smoothened by grav- 
ity. This temperature scale is related to the thickness 
of the superfluid-normal-fluid interface ~ 100 /^m by 
6Tg/^g^\dzTx\^ 1.273 MK/cni. 

To measure the critical singularity of the thermal re- 
sistivity or conductivity for heat currents smaller than 
100 nW/cm^ and temperatures closer than 15 nK to the 
superfluid transition, the experiment must be performed 
under microgravity conditions in space. Of course, for 
zero gravity where dzT\ = our theory presented in this 
paper works for all heat currents and never fails at any 
low Q. 



VI. THE TEMPERATURE PROFILE 

A. Numerical results 

Once the thermal conductivity At = At(AT, Q) has 
been determined, the temperature proflle T(z) is calcu- 
lated by solving the heat transport equation ( |5.l| ) as a 
differential equation which can be written in the form 



dT 

dz 



Q 



AT(AT,g) 



(6.1) 



where AT = T(z) - Tx{z). In Fig. | the resulting tem- 
perature profile of the su^erffuid-normal-ffuid interface is 
shown for Q = 1 /iW/cni . The several curves correspond 
to all three gravity conditions as indicated in the figure: 
vertical heat flow upwards and vertical heat flow down- 
wards on earth and heat flow in zero gravity in space. The 
temperature profiles T(z) are shown as solid lines. The 
spatially dependent critical temperatures Tx(z), which 
represent the gravity conditions, are shown as dashed 
lines. While in Sec. ^ we have found that At(AT, Q) 
does not depend on gravity for Q ^ 0.1 /iW/cm^, the 
temperature proflle T(z) is considerably influenced by 
gravity which is clearly seen in Fig. ^ because there are 
several solid lines. 

In Fig. U we have chosen the coordinates so that the 
superfluid-normal-fluid interface, where AT — T{z) — 
Tx{z) = 0, is located at z = and T(z) = To. For this 
reason all curves intersect with each other in the point 
(0,0). For z < the helium is normal fluid because 
AT = T(z) -Ta(z) > so that the sohd fine is above the 
respective dashed line. On the other hand for z > the 
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helium is superfluid because here AT — T{z) — T\{z) < 
so that the respective sohd hue is below the respective 
dashed line. (Since the many curves in Fig. ^ may be 
confusing, one should have in mind that only those solid 
and dashed lines should be compared with each other 
which belong to the same gravity condition.) Further- 
more, we have chosen a strongly enlarged temperature 
scale (in /iK) which resolves the very small temperature 
gradients in the superfluid region z > 0. For this reason 
in the normal-fluid region z < the solid line in Fig. |^ 
(temperature profile T{z) ) has a very steep gradient and 
goes up nearly vertically for decreasing z. The shifted 
critical temperature AT\{Q), which here is —0.10 ^K, is 
located approximately at this point where the solid line 
has a "round corner" and the gradient changes from large 
negative to small negative. 
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FIG. 4. The temperature profile of the superfluid-nor- 
mal-fluid interface for the heat current Q = 1 /iW/cm^ for 
the three gravity conditions: vertical heat flow upwards and 
downwards on earth and heat flow in zero gravity in space. 
While the sohd lines correspond to the temperature profiles 
T{z), the dashed lines represent the spatially dependent crit- 
ical temperatures T\{z), which refiect the gravity conditions. 

In the normal- fluid region z < where AT > ATx{Q) 
the three solid lines fall all together into one line which 
indicates that here gravity is negligible. On the other 
hand in the superfluid region z > where AT < AT\{Q) 
all three solid lines differ from each other which means 
that here gravity is important. A general criterion can 
be found to distinguish the regions where gravity is im- 
portant and where it is negligible. In the differential 
equation ( |6.1[ ) gravity is included implicitly via AT = 
T(z) — T\{z) by the gradient of the critical temperature 
T\{z). Gravity will be important or not if the spatial 
dependence of T\{z) or of T{z) is dominating. We may 
define the Q-dependent temperature difference ATi(Q) 
by the condition |Vr| = |VTx|, which in terms of (S.l) 
can be written as 



At(ATi(Q),Q) 



dz 



1.273 /xK/cm . (6.2) 



For AT > ATi(Q) it is |VT| > |VTa| so that the 
heat current is dominating. On the other hand for 
AT < ATi(Q) it is I VT| < | VTa| so that gravity is dom- 
inating. Consequently, for temperatures AT sufficiently 
well above ATi(Q) it is |VT| < |VTa| so that gravity 
can be neglected while otherwise gravity effects are signif- 
icant. Thus, ATi{Q) may be viewed as the temperature 
which separates the heat current dominated region from 
the gravity dominated region. 

For the heat current Q — 1 /iW/cm^ in Fig. ^we find 
ATi{Q) = -0.14 fiK which is shghtly below but close 
to ATa(Q) = -0.10 ^K. Thus, in this case the inter- 
face between the superfluid and the normal-fluid region 
nearly coincides with the interface which separates the 
gravity dominated region from the heat current domi- 
nated region. This fact is clearly seen in Fig. § For 
other heat currents Q the situation may change, because 
ATi(Q) and AT\{Q) may be farther apart from each 
other. In Fig. | we plot ATi(Q) and ATx{Q) as func- 
tions of Q on a double logarithmic sca le. Our theoretical 
result for ATi(Q) obtained from ( |6.2| ) is shown as solid 
line. Furthermore, ATa(Q) obtained from (5.3) is shown 



as dashed line. For the heat currents Q > 0.1 /iW/cm^, 
for which our theory is valid in gravity, both ATa(Q) and 
ATi {Q) are negative. While in Fig. || the solid line repre- 
sents ATi{Q) for vertical heat flow downwards, ATi{Q) 
for vertical heat flow upwards will be slightly different. 
However, the difference is very small. It is smaller than 
the width of the solid line so that it can be neglected. 
Thus, the solid line in Fig. ^represents ATi{Q) for both 
heat flow directions with sufficient accuracy. 
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FIG. 5. The temperature shifts ATi(Q) (so lid l ine) and 
IS.T\(ff) (dashed fine), obtained from (3.2) and (5.3), respec- 
tively, as functions of the heat current Q in a double logarith- 
mic plot. ATi(Q), separates the gravity dominated region 
from the heat current dominated region, while ATa(Q) sepa- 
rates the superfluid region from the normal-fluid region. 
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For Q = 0.22 ^W/cm", in Fig. ^ the sohd hne and 
the dashed Hne intersect each other so that ATi{Q) = 
ATx{Q). For larger heat currents ATi{Q) is always be- 
low AT\{Q) so that the separation between the gravity 
dominated and the heat current dominated region is lo- 
cated always in the superfluid region or equivalently the 
superfluid-normal-fluid interface is located in the heat 
current dominated region. It turns out that for suffi- 
ciently large heat currents, say Q > 10 ^W/cm^, the 
superfluid-normal-fluid interface is nearly free from in- 
fluences of gravity, while gravity effects rise to a signifi- 
cant magnitude only far away from the interface in the 
superfluid region. However, since the interface thickness 
must be larger than the size of the thermometers, experi- 
ments to resolve the temperature profile of the superfluid- 
normal-fluid interface must be performed for very small 
heat currents Q < 0.1 /iW/cm^. For these small heat 
currents, on earth the interface would be strongly in- 
fluenced by gravity. For this reason, an experiment to 
measure the temperature profile of the interface is pre- 
pared-Lo be performed under microgravity conditions in 
spaceo 



B. Asymptotic formulas 

By using the asymptotic formulas of the thermal con- 
ductivity At (AT, Q) in Sec. 0.B, the differential equation 
(6.1) can be solved explicitly so that asymptotic formulas 
for the temperature profile T{z) are found in several re- 
gions. First of all we consider the heat flow in zero gravity 
where Tx{ z) = T\ is constant. In this case the differential 
equation ( |6.l| ) can be integrated easily by separation of 
the variables so that 



g-i AT(T'-rA,Q) dT' 



(6.3) 



In the normal fluid region f or T sufficiently well above 
T\ the asymptotic formula (5.4) may be inserted. Con- 
sequently, we obtain 



AT{z)^T{z)~Tx^K+{-Qz 



_n^U/[('i-2)i.] 



(6.4) 



for z <C with a certain constant K+. On the other hand, 
for T sufficiently well below T\ the asymptotic formula 
(5.14) may be inserted which implies 



Ar(z) = T{z) - Ta w -K^iQ^z 



3^U/[(3d-4)H 



(6.5) 



for z 3> with a certain constant /v_. Including the 
nonasymptotic effects of the dynamic RG theory, for the 
temperature profile T(z) we obtain the asymptotic for- 
mula 



T(z) 



Tx + K+{-Qz) 
Tx-K^{Q^z) 



1.79 
0.294 



for z < 
for z > 



(6.6) 



which is valid for the heat flow in zero gravity. This 
formula can be compared with the solid line in Fig. H 



for zero gravity. Clearly, the temperature profile T{z) 
decreases monotonically with increasing z. There is no 
lower bound for T(z) in the limit z — > cx). 

Asymptotic formulas for the temperature profile T(z) 
can be found also for vertical heat flows in gravity. The 
interface, which separates the heat current dominated 
region from the gravity dominated region, is located at 
a coordinate zi defined by AT(zi) = ATi{Q). For 
Q > 0.22 /iW/cm^ it is zi > so that the interface 
is located in the superfluid region. For z <C zi gravity 
is negligible. Consequently, the asymptotic formula ( |6.6[ ) 
remains valid for z <C zi where T\ = T\{z — 0) should 
be inserted. In the case of a vertical heat flow upwards 
another asymptotic formula can be found for z ^ zi. In 
Fig. 1^ it is clearly seen that in this case the temperature 
profile approaches a limiting value 



lim T(z)=Too(g) 



(6.7) 



which is below T\{z — 0). For z ^ zi the tempera- 
ture difference is approximately AT(z) « Tao{Q) —T\{z) 
so that its z dependence is governed by T\(z), i.e. by 
gravity. Then, by using the asymptotic formula At 
(— AT)^ '*(5^^ for the thermal conductivity we obtain the 
asymptotic formula 



T(z)«Too(Q) + A-ooQ'(z-zo)' 



for z > zi (6.8) 



where K^a > and zq < are certain constants and 
Too(Q) = Ta(z = 0) + zo9.Ta. 

On the other hand, for a vertical heat flow downwards 
the solid line in Fig. ^ approaches a straight line parallel 
to the dashed line for z zi so that the gradients d^T 
and dzT\ are nearly equal. The distance from criticality 
Ar(z) = r(z) — T\{z) approaches a constant value given 

by 



lim AT(z) = Ari(Q) 



(6.9) 



Consequently, for the temperature profile we find the 
asymptotic formula 

T{z)^Tx{z) + ATi{Q) forz>zi. (6.10) 



C. Dissipative region 

For t he v ertical heat flow upwards the asymptotic for- 
mulas ( |6.(^ ) and (3.8) indicate the existence of a dissipa- 
tive region which may be related to the dissipative r^ 
gion observed in the experiments by Liu and Abler, 
and by Murphy and Meyem. The superfluid transi- 
tion happens in two steps at z = and z = zi or 
equivalently at the temperatures AT — AT\{Q) and 
AT = ATi(Q). The region < z < zi, which corre- 
sponds to ATi(Q) < AT < ATx{Q) (in Fig. | the region 
between the solid and the dashed line) , may be identified 
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as the dissipative region, because here the teniperature 
profile T{z) has a finite gradient. On the other hand, 
the region z > zi which corresponds to AT ^ ATi{Q) 
(in Fig. I the region above the sohd fine) is the really 
supe rfluid region, because here the asymptotic formula 
(6.8) implies a nearly flat temperature profile with a 



nearly zero gradient. Thus, ATi{Q) may be interpreted 
as the transition temperature between the dissipative re- 
gion and the really superfiuid region, which should be 
related to the temperature shift 

ATe(g)=To,(g)-TA(z = 0) (6.11) 

measured in the experimentsEliji. We find that ATdQ) 
obtained from our theory is close to ATi{Q) and located 
slightly above the solid line in Fig. ^. From the slopes 
of the lines in the double-logarithmic plot we obtain ef- 
fective power laws Ari(g) and AT^Q) 

with nearly the same exponent x, which varies between 
X « 0.9 for 0.1 ^W/cm^ < Q < 1 A*W/cm^ and x w 1.25 
for Q > 10 AiW/cm^ 

While our theoretical predictions for the dissipati- 
ve region rSfitee qualitatively with the experimental ob- 
servationgj'Efcl, there are three major quantitative dis- 
agreements. First of all, the experimentally observed dis- 
sipation is much larger than the theoretically predicted, 
because for AT^Q) < AT < ATx(Q) the measured ther- 
mal resistivity — l/X^ oi Ref . || (crosses in Fig. ||) is 
much larger than the theoretically predicted (solid line 
in Fig. H). Secondly, the experimental ATdQ) ~ —Q^ 
does not agree with the theoretically predicted because 
the exponent Xexp = 0.81 of Duncan et al.El is consid- 
erably smaller than the theoretical exponent x > 0.9. 
Furthermore, our theory predicts a much larger spa- 
tial extent Az = zi jjif the dissipative region than it 
is in the experimentsEfQ. We find Az > 0.25 cm for 
Q > 10 /iW/cm^ which is about the sample size, while 
the experiments find a spatial extent Az much smaller 
than the sample size. 

The disagreements may possibly be due to the exper- 
iments of Refs—^,^^ because in the new experiment by 
Baddar et alO the dissipative region has not been ob- 
served in, t his f orm. We may use the experimental fit 
formulaO (^.15 ) for the thermal conductivity to c alcu- 
l ate t he related Ari,exp(Q) and Arc,cxp(Q) by ( [6^ ) and 
(6.11), respectively. As results we obtain power laws 
~ —Q^ with the exponent Xexp = 1-26. This exponent 
agrees quite well with our theoretical value x — 1.25 
for large heat currents. If we plot ATc.cxp(Q) in Fig- 
1^, the respective line would be parallel to the solid line 
for Q > 10 fiW/cm^ but located somewhat below the 
solid line. Thus, there remains a quantitative disagree- 
ment which is related to the quantitative disagreement 
of At. However, ATc_exp(Q) does also unt agree with 
the ATc{Q) of the previous experimentsEra. Thus, for 
the clarification of the disagreements also further exper- 
imental work is necessary. 

Since the thermal conductivity At (AT, Q) (solid line in 



Fig. 1^) is smooth and does not show any unusual behav- 
ior at ATi{Q) and ATc{Q), these temperatures are not 
properties of the helium. Rather, ATi(Q) and ATc{Q) 
are implied by gravity a nd o ccur when integrating the 
heat transport equation (5.1) or (3.1) to calculate the 
temperature nxpfile. Thus, we predict that in the experi- 
ment in spacecJ, where gravity is zero, the Arc(Q) of the 
experiments of Refs. ||, and || should not be observable 
and not be existent, while the AT\{Q) is expected to be 
found. 



D. Self organized critical state 

In gravity for vertical heat flows downwards the gradi- 
ents dzT and dzT\ are both negative. A situation may 
arise where both gradients are equal, dzT = dzT\, so 
that AT(z) — ATi{Q) is constant over a larger region 
in space. This state of the helium represents a self or- 
ganized critical (SQC) state which was considered theo- 
retically by OnokiEj and proposed for an^fixperiment by 
Ahlers and LiuLj. Recently, Moeur et al.El have realized 
the SOC state and measured the distance from criticality 
AT(z) = ATi{Q) as a function of the heat current Q for 
40 nW/cm < Q < 6 /iW/cm . The experimental re- 
sult agrees quite well with our theoretical prediction for 
ATi{Q), while our theory does not include any unknown 
adjustable parameters. This has been demonstrated in 
our previous rapid communicationtj for heat currents 
below 1.5 iiW/cm^ (see Fig. 3 therein). However, for 
Q > 1.5 fiW/cm deviations occur: for a given distance 
from criticality AT — ATi{Q) the related heat current 
Q is larger in the experiment than in our theory. This 
fact means that for Q ^4.5 /iW/cm^ the dissipation ob- 
served in the experimentJj is smaller than the dissipation 
predicted by our theory. 

The SOC state is stable in the following sense: the dis- 
tance from criticality AT(z) always converges to t he con- 
stant value Ari(Q) for large z according to ( |6.9| ). This 
fact is clearly seen in Fig. However, there are two pos- 
sibilities: AT(z) may converge to ATi{Q) either from 
above or from below. In the first case an interface be- 
tween normal-fiuid helium and the SOC state is found, 
while in the second case an interface between superfluid 
helium and the SOC state is found. In Fig. ^ the tem- 
perature profiles T(z) are shown for both cases, where 
the superfluid-SOC interface is represented by the lowest 
solid line. Both kinds of ir^erfaces were realized in the 
experiment by Moeur et al.Ll. 

Since AT(z) = ATi(Q) is constant, the SOC state is 
homogeneous in space so that it is an ideal system for 
theoretical and experimental investigations. In the Ap- 
pendix we evaluate the Green's function G(r, r') and the 
related quantities ng and Jg by assuming ri and Ar to be 
linear functions of the space coordinate r given by ( A12 ) 
and ( Al3| ). While in general this assumption implies an 
approximation, for the SOC state the assumption is ex- 
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actly satisfied, because ri and Ar are directly related to 
the temperature profiles AT — T{z) — T\{z) and Tx{z) 
which are constant and linear in z, respectively. 



VII. CORRELATION LENGTHS 

The RG theory includes a characteristic length defined 
by C = (mO"^! which is called the correlation length. 
Near criticality p«p-|theory yields the asymptotic result 
^ = ^oT'" whereElka v = 0.671 and = 1-45 x 10^^ cm. 
The correlation length depends on AT and Q indirectly 
via the RG flow parameter r — t{AT,Q) determined in 
Sec. 1^. In Fig. |g the correlation length ^ is shown log- 
arithmically as a function of AT for Q = 42.9 /xW/cm , 
i.e. the same heat current as in Fig. |. Clearly, ^ in- 
creases when AT approaches the superfluid transition. 
The nonzero heat current implies that f is a smooth func- 
tion of AT which has a maximum located at AT\{Q). In 
Fig. 1^ the maximum is clearly shown by the solid line, 
where its position at AT\{Q) is indicated by the arrow. 




-4 
AT [^K] 

FIG. 6. The correlation lengths as functions of AT for the 
heat current Q — 42.9 /iW/cm^. The correlation length ^ is 
shown as solid line. The dephasing length obtained from 
( [r.7| ) is shown as dashed line. The dotted line represents the 
characterist ic len gth scale ^2 of the temperature variations 
defined by (7.11). The arrow indicates ATa(Q). 



Correlation lengths can be observed in the equal-time 
Green's function G(r, r') because they represent the char- 
acteristic length scales for the decay of the Green's func- 
tion, when the separation of the two space points |r — r'| 
increases. In the normal fluid region for AT > AT\{Q) 
we find 



G(r,r')^0 for|r-r'|>f 



(7.1) 



which is valid for zero and nonzero heat currents Q. In 
the superfluid region for AT < AT\{Q) the decay hap- 
pens in two steps according to 



G(r, 



^2 gjk(r- 





for C < |r- r'l < 
for |r-r'| > 



(7.2) 



where is a second correlation length which will be de- 
fined below. In this case the decay is strongly infiuenced 
by the heat current via k and ^i. In thermal equilibrium 
where Q = and k = 0, the long-range order implies 
a nonzero order parameter (ip) = rj so that f i is infinite 
and G(r, r') does not decay to zero but only to the ab- 
solute square of the order parameter rj^ . In the previous 
theory of Ref . ^ the heat flow is assumed to be metastable 
where dissipation by creation of vortices was neglected. 
Thus, in this theory long-range order is preserved also for 
nonzero Q so that = rye*''"" is again nonzero, is 

infinite, and the Green's function does not decay to zero. 

However, in the present theory vortices are included 
implicitly so that {ip{r)) = 0, is finite, and G(r, r') 
decays to zero eventually for large distances. The second 
correlation length represents a dephasing length of the 
order-parameter field tlj{r) which is equal to the average 
distance between the vortices. An approximate formula 
for_£i can be extracted from the integral representation 
( A28[ ) of the Green's function in the Appendix. In the 



superfiuid region for AT sufficiently below AT\{Q) the 
exponential factor exp{— afi — a^s} as a function of a 
has a sharp maximum located at 



«o = (-C)^/V[3^/^(-fi)] 



(7.3) 



where ri and C 



-1/3 



fi are negative. Thus, in the 
second line of (A2S) we may replace a approximately by 
ao so that the integral can be evaluated by (A29). As a 
result we obtain the Green's function 



G(r,r') 
where 



exp{zk(r-r')-(r-r')V2C?} , (7-4) 



k = aa{F/8-fw')h 
is the wave vector and 

a = (2«o)^/' 



(7.5) 



(7.6) 



is the dephasing length. Clearly, Eq. (7.4) implies the 
asym ptoti c de cay formula (7.2) where 77^ — rig. Inserting 
(7.3) into (7.5) we find that k is exactly the average wave 
vector ( A41 ) which we have calculated in the Appendix. 



By using ri = ^ ^pi and pi = cr-'/^C we find that the 



dephasing length (7^) can be written in the form 



6 



1/4-1/6 



(7.7) 



This formula is suitable for a numerical evaluation of ^1 
because C and a are two of the seven dimensionless pa- 
rameters, which are determined in Sec. IV. In Fig. |^ 
the dephasing length is shown as dashed line. We 
note that ^1 is defined only in the superfiuid region for 
AT < ATx{Q). It turns out that ^1 is always larger 
than ^. While the correlation length ^ decreases with 
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increasing distance from criticality, the dephasing length 
^1 increases. 

By using the asymptotic formulas of Sec. ^.B an 
asymptotic formula for £i c an be deri ved. Eliminating 
the amplitude A from (5.1C) and (5.11) we obtain 



r-l/6 



(-0 



1/2 



8u[t] QS^d' 



Consequently, from (7.7) we obtain 

1 1/4 Ad goksTx 



8u[t] Qe 



d-2 



(7.8) 



(7.9) 



Since C is nearly constant, the leading dependences on 
AT and Q are governed by 



6 - Q-'^ 



-l^-(d-2) _ ^-1 



Q-\-AT) 



{d-2)u 



(7.10) 



This asymptotic formula clearly indicates that the de- 
phasing length £i diverges in the limit Q ^ 0. Thus, in 
thermal equilibrium where Q = Q the second correlation 
length £i is infinite as expected, which means that in this 
case vortices are not present. 

In the Appendix we have calculated the Green's func- 
tion G(r, r') and the related quantities Ug and Jg approx- 
imately by linearizing the parameters ri and Ar l ocally 
with respect to the space variable r according to (A12) 
and (A13). Since ri and Ar are related to the temper- 
ature profiles Ar(z) = T{z) - Tx{z) and T(z), respec- 
tively, and since T\{z) is always a linear function in z, the 
validity of the approximation is proven if the curvature 
&lT of the temperature profile is negligible compared to 
the gradient dzT. For this purpose we define 



^,^\d.T\/\dlT\ 



(7.11) 



which may be viewed as the characteristic length scale 
for the variations of the temperature profile T{z) with 
respect to the space coordinate z. In Fig. ^ the charac- 
teristic length ^2 is shown as dotted line. The curvature 
of T{z) is negligible if £2 is considerably larger than the 
intrinsic correlation lengths £ and £1 of the helium. In 
Fig. 1^ the dotted line is considerably above the solid line 
and the dashed line for temperatures AT sufficiently far 
away from lS.Tx{Q). Consequently, the criterion for the 
validity of our approximation is satisfied in the normal- 
fluid region and in the superfluid region sufficiently far 
away from the superfluid-normal-fluid interface. Close 
to the interface all three correlation lengths £, £1, and 
£2 have the same order of magnitude so that here the 
curvature of T{z) is important and our approximation is 
strictly speaking invalid. However, since nothing unusual 
is observed here, we believe that our theory yields reason- 
able and reliable interpolations for the physical quantities 
close to the interface. 

The dotted line in Fig. |^ represents the characteristic 
length £2 for a heat flow in zero gravity. Similar results 
for £2 are obtained if we insert the temperature profile of 



a vertical heat flow in gravity. Deviations are expected 
for temperatures AT % ATi{Q). In gravity we find a 
somewhat smaller £2 , except for the downwards heat flow 
if AT is close to ATi (Q) . In the latter case £2 diverges in 
the hmit AT ATi(Q) which means that for the SOC 
state the approximation is exactly valid. Eventually, it 
turns out that the criterion for the validity of our approx- 
imation for calculating G{r, r'), Us, and Jg is not affected 
significantly by gravity. 



VIII. ENTROPY AND SPECIFIC HEAT 

In model F the field variable m(r, t) represents a fluc- 
tuating entropy density divided by fee- For this reason 
the local entropy density is defined by 



S = So + kB (m) 



(8.1) 



where 6*0 is a constant. Thus, the average ( to) m ust 
be evaluated. From the free energy functional (2.3) we 
derive 



(-) 
\5m/ 



-7o(|V'P)-/io 



We insert this quantity into (3.29), replace (IV'I ) 



and obtain 

ro{z) = To(z) + 2xo7o'^s + 270 (to) 



•2) 



(8.3) 



From ( |3.3lD we obtain Us — (4uo)^^[ri(z) — 7*0(2:)]. We 
resolve (|8.3D with respect to (to) and then obtain the 
entropy 



S = S< 



'o + A:b(27o) ^{-ro(z) + ro(z) 



+ ^i^oiz)-niz)]] . (8.4) 

Within our approximation in Sees. || and |l| we find 
ro{z) = ri(z) = at the critical point (AT, Q,^) = 
(0,0,0). Hence 



5a = 5*0 - fee (270) Vo(z) 



(8.5) 



is the entropy density at T = Ta in thermal equilibrium 
and zero gravity. We define the entropy difference AS* = 
S — S\ as the deviation fr om the entropy at criticality 
5a. Then, from 



and 



A5 



^o(^ ) 
' 27 



-{' 

I- 



Xo7o 
2uo 



5) we obtain 
ri{z) 



ro{z) 



(8.6) 



While 5a is a constant, A5 = AS {AT, Q) is strongly in- 
fluenced by the critical fluctuations near the superfluid 
transition. For this reason, we must renormalize the en- 
tropy A5 and apply the RG theory. In analogy to-the 
field variable TO(r, t) the entropy is renormalized byllj 



AS = (xo^™)'/'A5r, 



•7) 
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In (S.6) we replace the bare model-f par a mete rs by 
t he renormalized counterparts by using (4.1), ( [4.5| ), and 
(4.6). As a result we obtain the renormalized entropy 



27 



2_ 

2u 



1 



r{z) 



Within the Hartree approximation all Z factors cancel, 
where an expansion with respect to the coupling param- 



eters is not necessary. By resolving (4.25) with resp ect 
to Xo^m we obtain the renormalization factor of ( ^.7| ) as 

{XoZ„,f'^ - (Arf/)i/2 [2r7]-i . (8.9) 



by the di mens ionless variables 
We apply the 



We replace r(z) and ri(z] 
p and pi accor ding t o (4.11 ) and (4.13). 
RG theory to (^.7|)-(pJD by replacing ^ 
u — > u[t\, 7 — > 7[t], etc., so that now all dimensionless 
coupling parameters depend on the RG flow parameter 
T. Combining the resulting three equations together we 
eventually obtain the entropy difference 



AS- 



1 



1 



■7m- 



2u[t] 



1-^ 



} . (8.10) 



Eq. 



.IC ) is the final formula which can be used for nu- 



merical evaluation of the entropy density S 
There are no adjustable parameters present in this for- 
mula. The needed dimensionless parameters r, p, and pi 
and the correlation length 1^ — ^(r) were determined in 
Sec. [V, We note that Eq. ( ^.10[ ) was derived within the 
Hartree approximation. It is not restricted to the phys- 
ical situation considered in this paper where the helium 
is influenced by a heat current and gravity. The formula 
can be applied also to other physical situations as e.g. 
rotating helium whenever the Hartree approximation is 
used. 

The entropy AS* is a static quantity so that nonasymp- 
totic effects of the dynamic RG theory are v ery small. For 
this reason, close to criticality Eq. ( 8.1C| ) can be eval- 
uated asymptotically by using the asymptotic formulas 
^ = S^qt'", u[t] w m* = 0.0362, and 



4i/^ 



7 m 



(1 - 



.11) 



where v = 0.671 and a = 2 — dv = —0.013 (see Refs. 
|l| and |3^). Eq. < ^T^ is obtained by solving the RG 
equation— for 7[t] asymptotically. Since 7[t] is a known 
functiorO, 6 is a known constant (which actually is close 
t o un ity). Now, inserting the asymptotic formulas into 
10) we eventually obtain the entropy density 



S^Sx + t{B + A[ {^vja) + E\u*\ ]r-"} (8.12) 
where 



t = Tp = ^T/Tx , 
E[u\ = {2u)-\l-pi/p] 

B = A{-Av/a)h . 



(8.13) 
(8.14) 
(8.15) 
(8.16) 



The specific heat C is obtained by differentiation of 
the entropy S with respect to temperature according to 

dS__dS_ 



C^Tx 



.17) 



From our experience the best way to calculate the specific 
heat is first to calculate the entropy S by ( B.1C| ) or ( ^.12 ) 
and then to deterrnine the specific heat C by numerical 
differentiation via ( S.17 ). (Th e alternative way, to differ- 



entiate the bare entropy (8.4) first with respect to ro(z 



and then to apply the RG theory, is less reliable and 
yields artifacts, so that it should not be used.) While 
the constants A and B are given by ( p. 15 ) and ( ^.16 ), 
we can alternatively determine these constants by fitting 
the specific heat in thermal equilibrium to the data of 
the newest experimentES, which was performed in micro- 
gravity in space. In this way we obtain A — 2.22 J/molK 
and B — 456 J/molK. Since in the experiments the mo- 
lar specific heat is measured, we have multiplied the con- 
stants by the molar volume Vx , which for saturated vapor 
pressure isEJ Vx — 27.38 cm'^/mol. 
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FIG. 7. The specific heat Cq{/^T, Q) as a function of AT 
for the constant heat curre nt Q = 42.9 uW/cm^. Our theo- 



retical result obtained from (3.12) and (S.17) is shown as soM 
line. The dashed line represents Cq of the previous theoryEj, 
where vortices were neglected. For comparison we have plot- 
ted the specific heat in thermal equilibrium for Q = as 
dotted line. The arrow indicates AT\{Q). 

The specific heat C depends on the thermodynamic 
variable which is kept constant when perfor ming the dif- 
ferentiation with respect to temperature in (B.17). Since 
we consider liquid "^He in the presence of a constant heat 
flow, the heat current Q is the natural variable which 
should be kept constant. For this reason we calculate 
Cq at constant Q. In Fig. Cq = Cq{AT, Q) is plotted 
as a function of AT for Q = 42.9 /iW/cm^, i.e. the same 
heat current as in Figs. || and ||. Our theoretical result 
is shown as solid line. For comparison, the specific heat 
in thermal equilibrium at Q = is shown as dotted line. 
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While in thermal equilibrium the specific heat is singular 
at AT = 0, for nonzero Q we find a smooth curve for Cq 
which exhibits a strong maximum located at ATx((5). 
We note that the temperature T{z) is space dependent, 
so that the specific heat Cq{z) is also space dependent 
and must be interpreted as a local quantity, j.. 

Within the framework of the previous theoryo, the spe- 
cific heats Cy^ a nd Cq = Cj^ were calculated and a_aimi- 
lar formula like ( ^.12| ) was olatained for the entropj^Zl. In 
Fig. 0the specific heat Cq of the previous theory is shown 
as dashed line. Since the heat fiow is metastable below 
and unstable above AT\((5), the dashed line is defined 



I heat Cq is en- 
lat ATA(g). 



only for AT < ATx{Q). The specif 
hanced by the nonzero Q and divergest 
periments to measure Cq at constant Q are in progresg23. 
An enhancement of the specific heat by a nonzero heat 
current was found just recently. 
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FIG. 8. The specific heat Cvt(AT, Q) as a function 
of AT for the self organized critical (SOC) state where 
the temperature gradient is fixed by gravity according to 
dzT — dzT\ = —1.273 ^K/cm. Th e sol id line repr esents 



our theoretical result obtained from (S.12) and (8.17). For 
comparison the specific heat in thermal equilibrium is shown 
as dotted line. 

The self organized critical (SOC) state represents an 
ideal system for measuring thermodynamic quantities 
like the specific heat because it is homogeneous in space 
over a large region. In this case the temperature gra- 
dient is fixed by gravity according to dzT = dzTx = 
— 1.273 fiK/cm. Thus, in the SOC state the specific heat 
CvT at constant tempera ture gradient will be measured. 
We calculate Cvt from 



17) by inserting the entropy 



12) and keeping dzT — dzT\ constant when perform- 



ing the numerical differentiation, while the heat current 
Q is varied appropriately. The result is shown in Fig. ^ as 
solid line. For comparison, the equilibrium specific heat 
is shown as dotted line. Clearly, the nonzero temper- 
ature gradient dzT = dzT\ of the SOC state induced 
by gravity implies a rounding of the critical singular- 



ity. The solid line is smooth and exhibits a maximum 
at ATg = -45 nK. While in Fig. ^ the maximum of Cq 
is very strong and enhanced, here in Fig. ^the maximum 
of Cvt is moderate and just represents a smooth and 
round interpolation of the equilibrium specific heat near 
criticality. The temperature scale in Fig. || indicates that 
nano-Kelvin resolution is sufficient for a measurement of 
Cvt- Thus, the rounded temperature dependence of 
Cvt should be accessible by present-day experiments. 
Since gravity is needed for the realization of the SOC 
state, the experiment should be performed on earth. 



IX. VORTICES AND MUTUAL FRICTION 

While our theory does not include vortices explicitly, 
here we present some arguments, which support that our 
theory includes the effect of vortices implicitly. In Sec. 
[v] we have calculated a finite thermal conductivity At 
for AT < ATx{Q) which implies dissipation of the heat 
current Q in the superfluid state. Since the heat is trans- 
ported convectively by the superfiuid-normal-fiuid coun- 
terfiow, a superfluid current with a velocity Vg is induced 
in the opposite direction of the heat flow. The superfluid 
velocity Vg is related to the phase of the order parame- 
ter field tp{Y,t), so that the dissipation of the superfluid 
current is necessarily related to the creation of vortices. 
Thus, since our theory describes dissipation, it must in- 
clude vortices in somCr-svay. 

Gorter and MellinkliJ investigated mutual friction of 
the counterflow experimentally and proposed the mutual- 
friction force density 



(9.1) 



where A is the so called Gorter-Mellink coefficient which 
may be a function of temperature but which should be 
independent of the velocities Vs and v-o_. The force den- 
sity ( |9.l[ ) was added to the hydrodynamic equations of 
the two-fluid model (see (8) and (9) in Ref. |l5|). For a 
stationary counterflow the relation 



PnPs 



s\ dzT = Ap^ps{vs 



(9.2) 



was found, where the pressure gradient is neglected close 
to criticality. This relation implies dzT ~ [vs — Wn)^ 
and means that the temperature gradient induced by 
the mutual friction is proportional to the third power 
of the counterflow velocity. -By experimental and the- 
oretical considerations Vinenllj showed that the ansatz 



(9.1) is related to a turbulent superfluid flow which gen- 
erates a tangle of vortices which imply the mutual fric- 
tion. A statistical theory for the superfluid turbulenceuin 
a homogeneous counterflow was develo ped by Schwardl3. 
This latter theory supports (9.1) an(^ |9.2[ ). The Gorter- 
Mellink coefficient A was calculatedlj for temperatures 
in the interval 1.2 K < T < 2.05 K and agreement with 
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the experimentsESEJ was found (see also the review by 
ToughEl). 

Now, here we show that the ansatz of Gorter and 
Mellink (|]l]) can be derived from mod el F by our ap- 
proximation. To do this, we resolve (3.2) with respect to 
A, obtain 



A 



(9.3) 



and insert the results of our calculations. The coeffi- 
cient A must be proven to be independent of Q and 
Vs — v-o_. Near criticality t;n can be neglected because 
it is much smaller than Vg. The superfluid velocity Vs is 
related to the wave number k of the order parameter by 
Vs = hk/mi. The entropy|-Ber mass sx is related to the 
model-F parameter go bju s\ = {?i/m4){go/T\). The 
temperature gradient dzT is related to t he thermal con- 
ductivity At by (O). Then, from ( |9.3|) we obtain the 
Gorter-Mellink coefhcient 



A 



/m4\'^ go Q 



(9.4) 



\ h J pxTx Axfc^ 

In the previous theoryS the heat current Q was calculated 
as a function of the wave number k where the heat flow 
is metastable and dissipation by creation of vortices is 
neglected. In the superfluid region sufficiently well below 
AT\{Q) the dissipation is small so that the result for 
Q — Q{AT, k) of Rcf. H] may be used to replace k by Q. 
In this region, Q — Q{AT,k) is approximately a linear 
function of k given by 



Consequently, we obtain 



k 



(9.5) 



A 



9o 



V h J pxTx AtQ2 



(9.6) 



The last factor {• • •} does not depend explicitly on Q. It 
only depends on the RG flow parameter r via ^ = ^o'''"'^ 
where u\t\ k, u* = 0.0362. Now, inserting the asymptotic 



formula (5.13) for the thermal conductivity At we clearly 
see that the heat current Q cancels. Neglecting the one- 
loop contribution in the heat-current formula ( |9.5[ ) we 
obtain the Gorter-Mellink coefficient 



A 



go vT2 47 [T]w'[r] 

^n)px (-c)3/2 f[t] 



re 



(9.7) 



which does not depend explicitly on Q. A weak indirect 
Q dependence is included in the dimensionless variable 
C, which represents the nonasymptotic effects of the dy- 
namic RG theory. In leading order we find the asymp- 
totic formula 



(9.8) 



This result proves that within the approximation of our 
theory the ansatz of Gorter and Mellink is correct and 
can be derived from model F. 



Eq. (9.6) can be used for an explicit calculation of 
the Gorter-Mellink coefficient A = A{AT, Q) as a func- 
tion of AT = T — T\ and Q, where our result for 
At = Xt{AT,Q) from Sec. |^ is inserted. While most 
constants and parameters are known, we additionally 
need m/^/Ti — 6320 s/cm^ and the density p\ of the he- 
liujBi at T\, which at saturated vapor pressure is given 
byO px = N^^m^/Vx = 0.146 g/cuY. In Fig. | our the- 
oretical prediction for A is shown versus AT double log- 
arithmically for several heat currents Q which change 
by a factor of 10 between each curve. The lowest and 
most left curve corresponds to the smallest heat current 
Q = 1 pW/cm^, while the highest and most right curve 
corresponds to the largest heat current Q = 1 mW/cm . 
Since the helium is superfluid only for AT < ATx{Q), all 
curves have an endpoint on the left-hand side. 
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FIG. 9. The Gorter-Mellink coefficient A obtained from 
( |9.(| ) as a function of temperature for several heat currents 
between Q — 1 pW/cm^ and 1 mW/cm^ (curves from left to 
right, heat currents increase by a factor of 10 respectively). 
The left ends of the curves are located close to the superfluid 
transition at ATx{Q). 

If the cancellation of Q in (^^) is perfect, then the 
curves in Fig. ^ would lie all on the same line. This, 
however, is not the case. The nonasymptotic effects of 
the dynamic RG theory imply that A{ AT, Q) depends 
weakly on Q. While the heat current Q is varied over 9 
decades, A changes by a factor of 3 to 10 or by 0.5 to 1.0 
decades. Consequently, we approximately find 



(9.9) 



with an exponent y between 0.05 and 0.1^ This result 



slightly modifies the mutual friction force (9J^) of Gorter 
and Mellink into 
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(9.10) 



with the exponent (j) — 3 + y between 3.05 and 3.1. 

The Gorter-MeUink coefficient A-was measured as a 
function of temperature by Vinenli3. The value A « 
200 cms/g was obtained for AT « —0.1 K, which can 

be extrapolated to A ^ 600 cms/g for AT 0.01 K. 

Similar values for A were obtained also in later experi- 
ments, which are reviewed in Ref. |42[ On the other hand, 
the lowest values, that our theory predicts, are about 
A - 10^ cms/g for AT 0.01 K. Thus, our theoret- 
ical prediction for A is about one or two decades larger 
than the experimental values. Furthermore, A can be 
extracted from the experimental da ta of Baddar et al.EII 
by inserting the power-law formula ( 5.15| ) for the thermal 
conductivity At into (9.6). Again, our theoretical predic- 
tion is about a factor of 20 larger than the experimentally 
observed values. Thus, we conclude that our theory has 
a tendency to overestimate the magnitude of the dissipa- 
tion and mutual friction due to creation of vortices. On 
the other hand, the experimental data of Ref. |l^ impl y 
the exponent cj) = 3.53 for the mutual friction force ( 9.10 ), 
which is considerably laxger than the exponent proposed 
by Gorter and Mellinkll3 and obtained from our theory. 

The discrepancies are possibly due the number of vor- 
tices in the helium in the presence of a homogeneous heat 
flow Q, because the number of vortices appears to be an 
unknown and uncontrolled quantity in theory and exper- 
iment. In rotating helium mutual friction can by stud- 
ied in a much more controlled fashionH, because in this 
case the number of vortices is related to the rotation fre- 
quency. Here the mutual friction is described by the Vi- 
nen coefficientsc3 B and B' . These coefficients have been 
measured experimentally (see Ref. H for review) and 
also calculated within model F in renormalized mean- 
field theoryc3 by considering the motion of a single vortex 
line. We have applied our theory (the Hartree approx- 
imation combined with the RGL-theory) also to rotating 
helium and find the coefficientala B — (4m4/?i)r'[r] and 
2 - B' = (4TO4/^)r"[r]. While these results are indeed 
simple, they have the same order of magnitude than the 
coefficients B a,nd,B' of the previous theoryc3 and also of 
the experimental^. The agreement is better for B than 
for B' . Thus, we conclude that our present theory indeed 
includes the effects of vortices. However, since a simple 
approximation is applied, the Hartree approximation (see 
Fig. 1^), discrepancies are expected. Nevertheless, for the 
Vinen coefficients B and B' the discrepancies are much 
smaller than for the Gorter-MeUink coefficient A. 



as the Hartree approximation, combined with the RG 
theory. While above T\ the previous theory of Ref. || 
is recovered for a heat flow in normal fluid ^He, below 
T\ in the superfluid state our theory predicts dissipation 
of the heat current and mutual friction of the related 
superfiuid-normal-fluid counterflow. _ 

We derived the ansatz of Gorter and MellinkEj for the 
mutual friction force and found several indications that 
our approach includes vortices indirectly. However, our 
approach appears to overestimate the magnitude of the 
dissipation by vortex creation considerably compared to 
the experimental observations. This discrepancy is prob- 
ably due to the number of vortices in the superfluid he- 
lium which appears to be an uncertain and very sensitive 
quantity strongly influenced by the kind of the approxi- 
mation and also by the experimental conditions. Further 
theoretical and experimental work is necessary to clarify 
the discrepancies. 

Besides the correlation length ^ = ^qt"", which is the 
conventional length scale of the RG theory, in the su- 
perfluid state we flnd a second characteristic length 
which describes the decay of the correlations by dephas- 
ing of the order parameter field ■0(1") caused by vortices. 
The dephasing length may be viewed as the average 
distance between the vortices. It is much larger than ^ 
(see Fig. ^) and approximately given by « aQ^^r^" 
where r = -2AT/Tx and a « 1.0 mW/cm. Thus, for 
superfluid ^He in confined geometries we expect novel ef- 
fects when the dephasing length is as large as the ge- 
ometry and hits the boundary walls. Future theoretical 
and experimental investigations should study the influ- 
ence of vortices in superfluid ^He on flnite-size effects and 
on boundary effects. Unexpected results for the thermal 
conductivity in confined geometries and for the Kapitza 
resistance may possibly be found, which are caused by 
the second characteristic length ^i. 



ACKNOWLEDGMENTS 

I would like to thank Professors G. Ahlers, V. Dohm, 
R. Duncan, and H. Meyer for valuable discussions. Fur- 
thermore I would like to thank Drs. M. Lee, U. Israelsson, 
and D. Strayer for the invitation to the 1998 NASA/JPL 
workshop on fundamental physics in microgravity, held 
in Oxnard, CA, where many question about the relation 
between theory and experiment were discussed and clar- 
ified. 



X. CONCLUSIONS 



APPENDIX: INTEGRAL REPRESENTATION 
FOR THE GREEN'S FUNCTION 



We have presented a renormalization-group (RG) the- 
ory based on model F for liquid "^He near the superfluid 
transition in the presence of a heat current Q and gravity. 
The fundamental concept is a self-consistent approxima- 
tion, which in quantum many-particle theory is known 



In Sec. Ill and in Ref. || we have evaluated Us — (IV'P) 
and Jg = (Im[^/'* V"!/)]) explicitly, whi ch ev entual ly ar e 
expressed in terms of an integral by ( 4.27 ) and ( 4.28 ). 
While for our purposes we only need rig and Jg, it is also 
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possible to evaluate the complete matrix Green's fu nctio n 
G explicitly. This Green's function is defined by (3.18). 
First of all, we need the inverse of the operator matrix 
(3.17) which is given by 



Since we intend to evaluate the renormalized Green's 



function, we use (4.4) for the operator L, where all model- 
F parameters are replaced by the renormalized parame- 
ters. 

Here we consider the Green's function (ipip*) which 
is given by the lower right element G22 of the matrix 
Green's function G. From (3.18|) and (Al) we obtain 



(^(r,i)r(r',0> =4r'L-i(L+)-i<5(r-r')<5(t-0 . 

(A2) 

We represent the inverse operators L^^ and (L+)^^ as in- 
tegrals of exponential functions so that the Green's func- 
tion is rewritten as 



(V'(r,i)V'*(r',t')) =4r' / da j (1/3 e""^ e'^^^^ 



Jo 
X S{r - r') d{t - t') 



(A3) 



We decompose the operators in the exponentials into 
time-dependent and space-dependent parts according to 



— aL 
-pL+ 



-adt 



aA , 
(3B , 



where 



A = -{T[n-V^]-z{g/2^)Ar} , 
B = -{T*[n-V^]+i{g/2^)Ar} 



(A4) 
(A5) 



(A6) 
(A7) 



Since dt commutes with the space-dependent operators 
A and B, we obtain 



(V'(r,i)V'*(r',t')> ==4r' / da d/3 
Jo Jo 



xe 



aA 3B 



<5(r-r') 

xe(-"+^)^'(5(t-t') 



(A8) 



The time-dependent factor is evaluated by the Taylor se- 
ries according to 



,(-«+/3)9t s{t - t') = 5{-a + (3 + t-t') 



(A9) 



The space-dependent factor is more complicated and can 
be evaluated by using the formula 

e^e^ = exp{A + B+ 

+ [A,B]]-[B,[A,B]]) 
-2.[B,[A,[A,B]]] + ---} . (AlO) 



1. The equal time Green's function 

For simplicity we set t' = t, and from now on we con- 
sider only the equal time Green's function. The inte gra l 
over /3 can be evaluated easily by the delta function ( A£ ) 
so that we obtain 



G(r,r')-(^(r,t)V'*(r',i)) 

= AT' / da e"^ e"-^ 5{y - r') . (All) 
Jo 

For a successful calculation of the space-dependent inte- 
grandthe series in the exponential on the right-hand side 
of (IaTcI ) must be finite, which means that only a finite 
number of the commutators may be nonzero. For this 
reason, we assume ri and Ar to be linear functions of 
the space coordinate r given by 



n = ai 
Ar = a H 



f bir , 
br , 



(A12) 
(A13) 



where oi and a are constants and bi and b are constant 
vectors. In general, ri and Ar are nonlinear functions of 
r. In this case ri a nd A r must be linearized locally so 
that Eqs. ([aT^ ) and (|A13|) are taken as an approximation 
where bi and b are the respective gradients at the space 
point R = ^(r -f r'). In the main text of this paper we 
assume ri(z) and Ar{z) to be functions of the coordinate 
z only so that the gradients bi = r[ez and b = Ar'e^ are 
vectors in z direction. However, this latter assumption is 
not necessary here. 

Now, Eqs. (A12) and ( A13| ) imply the commutators 
[n, V^] = -2biV and [Ar, V^] = -2bV. Consequent- 
ly, we find 

[A,B]^-t{gr'/^)2hV , (A14) 
[A,[A,B]]^ -2i(gr77){rbib-z(g/27)b2} , (A15) 
[B, [A, B] ] = -2z(gr77){r*bib + iig/2j)h^} . (A16) 



The double commutators ( A15| ) and ( A16) are c numbers 
because bi and b are assumed to be constant. Hence, all 
higher-order commutators are zero so tha t the series in 
the exponential on the right-hand side of (A10|) is f inite. 
Event ually, inserting the commutators (A14)-(A16) into 
(AlO) we obtain 



aA aB 

e e 



= exp|-2r'a[ri - V^] - i (2r'a)^ 



47r' 



-bV 



-{2T'af [2— f ^bVi - 4f ^1 

12^ 'IT' J \A"fT' 



} ■ (A17) 



This formul a is e xact i f ri a nd Ar are linear functions in 
r given by ( [aT^ ) and ( [A1^ ). 

Next, we consider the operators C — — 2r'ri and D — 
2r'( V — zk)^ where k is a constant wave vector. We find 
the commutators 
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[C,D] = (2r')^ 2 bi(V-ik) , (A18) 
[C, [C, D] ] = (2r')^ 2 , (A19) 
[D,[C,D]]={) . (A20) 

Since the double commutators are c numbers, again the 
higher-order commutators are zero. Thus, from (AlO) we 
obtain the formula 



„aC „aD 



exp{-2r'a[ri - (V - ik)^] 



+ (2r'a)2bi(V - ik) + \{2T'afhl] 
We choose the wave vector 



k = V'a 



-h - ihi 



(A21) 
(A22) 



and then rewrite (A21) as 
goc gaZ) ^ exp(-2r'a[ri -V^]-i {2T'af-^hV 

-3^(2r'«f[b?+3(-^b)']}. (A23) 

Now, we clearly see that th e oper ators in th e exponentials 
on the right-hand sides of (A17) and (A23) are equ al up 
to an additive constant. Thus, from comparison of ( A17) 
and (A23) we obtain the relation 



aA aB 

e e 



where 



12 



F 



(A24) 



(A25) 



Here we have replaced F' , F" , and g by the dimensionless 
ratios w' = F'/A, w" = F'VA, and F = g/X. The first 
factor on the right-hand side of (A24) is just a constant, 
while the other e xpone ntial factors are operators. We 
insert the relation (A24) together with the oper ators C = 
-2F'(ai -t- bir) and D = 2F'(V - ik)^ into ( [aTi| ) and 
then obtain the Green's function 

G(r, r') = 4r' / da e-2(r'")^- g-2r'a(ai+bir) 

"'0 

xe2r'"(^-*)^5(r-r') . (A26) 

We substitute 2F'a a and evaluate the last two factors 
by using the Gaussian integral and the Taylor series: 



^Q(V-ik)" 



^(r 



j d'^u e-(4")"'"'-»ku j^j. _ r' + u) 

(A27) 



(47ra)-'^/2 
(47ra)-^/2 

(47ra)-'^/2 e»k(r-r')g-(4a)-i(: 



Thus, inserting the wave vector (A22) with the proper 



subsitution for a and with g/V replaced by F/w' , for 
the Green's function we eventually obtain the formula 



G(r,r') 



X exp< 



(47r)'^/2 ^0 
F 



da 

7^ 



exp{— ari — a'^ s\ 



{^a^Hr-r')-lir-r'f]. (A28) 



where fi = oi + ibi(r + r'). This formula is a simple in- 
tegral, there are no operators in the integrand any more. 
Clearly, the Green's function depends on the average co- 
ordinate R = i(r-l-r') implicitly via fi and explicitly on 
the relative coordinate Ar = r — r'. 



The integral (A28) has a similar but more general 
structure than the integral J-a{C,) in (4.28). Thus, it is 
possib le to obtain Ug and Jg from the Green's function 
([A28|). We find 



n3 = (|Vn-G(r,r) 
2 r°° da 







a' 



d/2 



exp{-Qri - a^s} (A29) 



and 



J. = (Im[^*V^]) = {2i)-\V - V] G(r, r') |^,^ 
F h da 



A-fw' (47r)'^/2 



;^d723T exp{-ari -as}. 

(A30) 



We identify bi 



e, and b = Ar'e, and find that s 



defined in (A2£) is closely related to a defined in (4.29). 
We substitute a — s^^/'^v, identify ( = s^^Z-^ri, and 
rewrite Ug and Jg as 



Js 



{AttY/ 
F Ar' 



A-fw' {Any/ 



,.('^-^)/«^(4-.)/2(C) 



(A31) 
(A32) 



where Jg = Jge^. By using the function ^a{X) defined 



in (4.27) and identifying 



1 



{AttY/- 



1 



F(l - d/2) 



(A33) 



which may be viewed as the definition of the factor Ad^ 
the formulas for and X can be rewritten as 



-l + £/2 



{X) r\ 



-e/2 



Js = 



F Ar' 1 
2^7 '2r/' e 



Arf(l-e/2) $,/2(X) r. 



-e/2 



(A34) 
(A35) 



where X = —s/rf and e — 4 — d. Finally, replacing 
the renormalized couplings by t he ba re mo del- F param- 
eters, we recover the formulas ( ^.24 ) and (3.25). Thus, 
we have derived the formulas of Sec. III.B for rig and Jg 



once again. We note that the fundam ental assu mptio n 
here is the linear form of ri and Ar in ( |A12| ) and (A13), 
which in general is an approximation. This a ssump tion 
implies the special structure of the integrals in ( A28) and 
in (Osl). 
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2. Fourier transformation of the Green's function 
and physical interpretation 

The natural space variables of the equal time Green's 
function G{t, r') are the average coordinate R — ^(r+r') 

and the relativ e coo rdinate Ar = r — r'. This fact is this we apply the formula (AlO) to the factor e°''^e^^ 



3. Concluding remarks 

We have evaluated the Green's function G22 = (V'V'*) 
only for equal times t' — t. However, it is possible to 
evaluate this Green's function also for t' ^ t. To do 



clearly seen in (A2S). Thus, we may perform a Fourier 
transformation with respect to Ar and define the Fourier- 
transformed Green's function G{Il, k) by 



G(r,r') = 



,ik(r-r') 



G{R, k) 



(A36) 



We apply this Fourier transformation to ( A28 ) and then 
obtain 



G'(R,k) = 2 / da exp{-afi - a^s} 



X exp-f— afk — a \ 



87W' 



(A37) 



where fi = Oi + biR. This Green's function is positive 
definite and may be viewed as a distribution function 
for the wave vector k. We may define the average wave 
vector (k)R at space point R by 



d'^k 

(2^ 



k G(R, k) 



d'^k 
(2^ 



G(R,k) . (A38) 



This formula can be rewritten in terms of the real-space 
Green's function G(r, r') as 

(k)K = {{2^)-'[V - V] G(r,r')}/G(r,r')L,^,^^ . 

(A39) 



Thus, by using ( |A29| )-( |A32| ) we obtain 

F b (-a^./2(C) 
-ri "~ 



(k)R = Js/?^s = -z 7 



8jw' 



l+c 



/2(C) 



(A40) 



so that Js = ris (k)R,. While Jg is finite for d < 4, rig 
is ultraviolet divergent for d > 2. As a consequence, for 
d = 3 dimensions the average wave vector (k)Fi is strictly 
speaking zero. In the superfluid state for temperatures 
well below T\ we may use the asymptotic formula (4.44) 
for the function J-a{(^). In this approximation the ultra- 
violet divergences are neglected so that we obtain a finite 
result for the average wave vector 



b (-0 



3/2 



F 



8jw' -n 31/2 



(A41) 



Since the wave vector k is related to the superflui d ve- 
locity Vs by Vg — Kk/m4, the Green's function ( A37 ) 
impli es a distribution function for the superfluid velocity 
and (A41) yields the average superfluid velocity (vs)r at 
space point R. 



in (AS). The eventual result is again an integral over a 
single variable a with, however, a somewhat more com- 
plicated integrand than (A28). Furthermore, the nondi- 
agonal Green's function G12 — {tpip) can be calculated. 
In this case there is only one inverse operator which 
implies only one integral over a. This integral is per- 
formed trivially by a delta function analogous to (^^^) 
so that eventually there will be no integral at all. The 
other nondiagonal Green's function G21 = G^2 is just the 
Hermitian conjugate. Finally, the upper left diagonal el- 
ement Gil is zero. Thus, we conclude that the complete 
matrix Green's function G defined by ( 3.18| ) can be eval- 
uated explicitly, supposed the parameters ri and A r are 
linea r func tions of the space variable r according to ( A12 ) 
and (|Al3|). 
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